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CHAPTER  I 


INTRODUCTION 

1.1  General  Aspects  of  the  Problem 

Over  the  past  two  decades,  convection  phenomena  induced  by 
body  forces  have  been  the  object  of  extensive  research  efforts.  A 
rather  broad  classification  into  two  types  of  problem  can  be  made: 
The  external  problem  which  considers  the  flow  and  heat  transfer 
over  a  body,  such  as  a  heated  rod  or  plate,  in  a  fluid  at  rest  and 
is  referred  to  as  free  convection.  The  internal  problem  which 
deals  with  the  flow  and  heat  transfer  in  confined  fluids  and  is 
referred  to  as  natural  convection. 

Natural  convection  in  completely  confined  fluids  occurs  in 
many  and  diverse  applications  and  has,  as  a  result,  been  receiving 
more  attention  recently.  The  study  of  natural  convection  flows 
that  are  completely  bounded  by  surfaces  leads  to  a  number  of  diffi¬ 
cult  problems.  Natural  convection  is  extremely  sensitive  to  the 
configuration  and  boundary  conditions.  Theoretical  analysis  is 
limited  because  of  the  inherent  coupling  and  nonlinearity  of  the 
basic  equations. 

For  large  Grashof  and/or  Rayleigh  numbers,  the  flow  will  ob¬ 
viously  be  of  a  boundary  layer  type.  For  external  problems  Prandtl 
boundary  layer  theory  yields  the  same  simplifications  that  are 
helpful  in  other  problems  and  the  region  exterior  to  the  boundary 
layer  can  be  determined  once  the  boundary  layer  flow  is  known. 


But  In  a  completely  enclosed  region,  the  region  exterior  to  the 
boundary  layer  will  actually  be  enclosed  by  the  boundary  layer  and 
will  form  what  is  called  the  "interior"  or  the  "core  region".  As 
the  core  region  is  completely  surrounded  by  the  boundary  layer, 
the  flow  in  the  core  is  not  readily  determined  from  the  external 
boundary  conditions  but  depends  on  the  action  of  the  boundary  layer. 

Hence  the  boundary  layer  and  core  are  closely  coupled  so  that  it  is 
very  difficult  to  predict  the  core  flow  pattern  a  priori.  Since 
the  analytical  boundary  layer  approach  does  require  a  priori  know¬ 
ledge  of  the  core  configuration  at  the  outset,  it  is  precisely  this 
coupling  that  constitutes  the  main  source  of  difficulty  in  obtain¬ 
ing  analytical  solutions  or  even  in  getting  qualitative  ideas  of 
the  flow  patterns  to  internal  problems. 

1.2  Review  of  Existing  Work 

Since  the  present  work  is  directed  towards  those  cases  in 
which  the  fluid  is  heated-from-side  and  is  mainly  concerned  with 
the  prediction  of  flow  patterns  in  the  core,  existing  work  will  be 
reviewed  in  detail  from  that  point  of  view. 

Most  of  the  earlier  works  were  concerned  with  natural  convec¬ 
tion  in  horizontal  cylinders  and  high  aspect  ratio  rectangles.  A 
comprehensive  review  of  earlier  existing  works  has  been  done  by 
Ostrach  [1].  More  recent  research  on  natural  convection  flows  in 
horizontal,  vertical  and  tilted  rectangular  enclosures  was  reviewed 
by  Catton  [2]. 

i 
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As  was  pointed  out  in  [1],  in  the  earlier  analytical  works, 
e.g.,  Ostrach  [3],  Batchelor  [4],  Weinbaum  [5]  and  Menold  and 
Ostrach  [6],  except  Gill  [7]  and  Hantman  [8],  the  core  configura¬ 
tion  was,  a  priori,  considered  as  isothermal  and  rotating  regard¬ 
less  of  the  thermal  boundary  conditions.  A  series  of  experiments 
by  Ostrach  and  Pneuli  [9],  Ostrach  and  Menold  [10],  Brooks  and 
Ostrach  [11]  and  Sbazebari  and  Ostrach  [12]  showed  that  the  de¬ 
tailed  nature  of  the  core  flow  depends  on  the  thermal  boundary 
conditions.  Experiments  applicable  to  the  heating-from-side  case, 
e.g..  Brooks  and  Ostrach  [11],  Sabzebari  and  Ostrach  [12],  Martini 
and  Churchill  [13],  Eckert  and  Carlson  [14],  Elder  [15],  etc. 
indicated  a  thermally  stratified  core  with  relatively  slow  flow. 
With  the  exception  of  the  work  of  Poots  [16],  all  the  numerical 
work,  e.g,  Wilkes  and  Churchill  [17],  Elder  [18],  De  Vahl  Davis 
[19],  MacGregor  and  Emery  [20],  Newell  and  Schmidt  [21],  Quon  [22], 
Taylor  and  Ijan  [23],  etc.  was  performed  with  the  experimental 
work  in  mind  and  was  in  essential  agreement  with  this  picture. 

The  failure  of  the  assumption  of  the  isothermal,  rotating  core  for 
heating-from-side  case  was  brought  out  especially  by  the  solution 
of  Menold  and  Ostrach  [6].  Although  a  rotating,  isothermal  core 
has  not  yet  been  experimentally  observed  for  heating-from-side 
case,  the  possibility  of  its  occurring  and  other  patterns  too  may 
not  be  dismissed. 

In  contrast  with  the  considerable  work  done  on  natural  con¬ 
vection  in  horizontal  circular  cylinders  and  high  aspect  ratio 
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rectangles.  It  is  rather  surprising  that  relatively  little  work 
has  been  done  on  natural  convection  in  enclosures  with  low  aspect 
ratio  even  though  this  problem  has  many  engineering  applications. 

Boyack  and  Kearney  [24]  first  obtained  numerical  results  for 
a  low  aspect  ratio  rectangle  and  emphasized  the  need  to  study 
natural  convection  processes  for  application  to  auxiliary  cooling 
equipment  of  high  temperature  gas  cooled  reactors  and  solar  collec¬ 
tors. 

Solan  and  Ostrach  [25]  discussed  the  need  for  understanding 
such  phenomena  in  order  to  improve  the  growth  of  crystals  by 
closed-tube  vapor  deposition. 

In  1973,  Klosse  and  Ullersma  [26]  considered  a  gas  in  a  two- 
dimensional  low  aspect  ratio  rectangular  enclosure  using  fully- 
developed  flow  as  the  base  and,  as  essentially  assumed,  an  S-shaped 
velocity  profile  was  found  in  the  central  region  of  the  enclosure 
with  vertical  flows  in  the  end  regions. 

Later,  Cormack,  Leal,  Imberger  and  Seinfeld  [27]  ^[29]  studied 
buoyancy  driven  convection  in  shallow  lakes  as  a  shallow  cavity 
problem  analytically,  numerically  and  experimentally,  respectively. 

In  the  analytical  work  of  Cormack  et.  al.  [27],  they  investi¬ 
gated  the  problem  when  the  aspect  ratio.  A,  goes  to  very  small 
value  with  arbitrary  but  fixed  Grashof  number,  Gr^,  by  the  matched 
asymptotic  expansion  method.  Their  results  show  that  when  A  -«■  0, 
the  flow  consists  of  two  distinct  regimes,  a  parallel  flow  in  the 
core  region  and  a  second,  non- pa rail  el  flow  near  the  end  of  the 
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cavity.  They  also  show  that  the  primary  driving  force  for  the  mo¬ 
tion  is  the  horizontal  temperature  gradient  in  the  core  and  the 
temperature  distribution  to  the  first  order  in  the  aspect  ratio.  A, 
in  the  cavity  is  strictly  linear  in  the  horizontal  direction.  Thus 
the  dominant  mode  of  heat  transfer  is  pure  conduction.  However, 
they  could  not  indicate  explicitly  the  upper  limit  of  Gr^  within 
which  their  solution  is  valid.  They  just  gave  an  approximate 
criterion  by  requiring  only  that  the  second  term  in  the  expansion 
for  Kj,  the  single  parameter  of  the  parallel  flow  solution,  be 
small  relative  to  the  first,  i.e.,  Gr^  Pr  A  £10  .  A  comparison 
of  numerical  solutions  of  Cormack,  Leal  and  Seinfeld  [28]  with  the 
asymptotic  theory  shows  excellent  agreement  among  the  analytical 
and  numerical  solutions  proved  that  A£  0.1  and  Gr^  Pr  A  £10  .  In 
addition,  the  numerical  solutions  also  demonstrate  the  transition 
from  the  parallel  flow  regime  to  the  boundary  layer  limit;  for  A 
fixed,  Gr^  -♦  <*>.  In  Imberger's  experiment  [29],  most  of  the  flow 

features  mentioned  in  numerical  work  were  also  qualitatively  pre- 

2  2  3 

sent  in  the  experimental  flows.  But  when  Gr^  Pr  A  becomes  nearly 
1011,  the  mid-depth  becomes  an  isotherm  and  there  is  a  slow  central 
circulation  throughout  the  whole  cavity,  which  is  absent  in  the 
much  larger  aspect  ratios.  Therefore  it  was  pointed  out  in  [29] 
that  it  was  uncertain  whether  the  limit  Ra^  -  is  different  for 
A  small  and  A  large.  In  this  respect,  it  was  also  pointed  out  that 
the  theory  developed  in  [27]  for  small  values  of  RaH  A  is  inappro¬ 
priate  in  the  sense  that  it  predicts  parallel  flow  in  the  core 
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region  for  all  values  of  the  Rayleigh  number. 

Bejan  and  Tien  [30]  considered  the  case  of  small,  fixed  A  with 
Ra^  variable.  They  developed  three  characteristic  flow  regimes 
with  respect  to  the  value  of  RaH  and  derived  a  set  of  analytical 
results  for  the  Nusselt  number  corresponding  to  each  flow  regime. 
They  assumed  the  existence  of  parallel  flows  in  the  core  in  all 
three  regimes,  but  this  is  not  true  under  certain  conditions. 

In  the  experiments  of  Ostrach  et.  al .  [31],  Al-Homoud  [32], 
and  Kamotani  et.  al .  [33],  many  different  core  configurations  were 
observed.  Therefore  their  predictions  are  of  limited  and  undefined 
validity. 

In  the  experiment  of  Ostrach  et.  al.  [31],  it  is  indicated 
that  for  the  range  of  aspect  ratios  and  Grashof  and  Prandtl  numbers 
covered,  i.e.,  0.05  £  A  ^  O.5.,  27.7  £  Gr^  £  106  and  0.72  £  Pr  £ 

3 

1.38  x  10  ,  the  flow  pattern  is  predominantly  unicellular  and  paral 
lei  for  the  aspect  ratios  of  0.1  and  smaller.  Secondary  cells  are 
observed  near  the  enclosure  ends  when  the  aspect  ratio  is  0.2  and 
(Gr^PrA)  is  of  the  order  of  10  for  Prandtl  numbers  of  the  order 

3 

of  10  ,  where  the  horizontal  boundaries  are  semi -conducting. 

In  Al-Homoud' s  experiment  [32]  with  water,  it  is  shown  that  in 

D  Q 

the  range  of  2  x  10  <  Ra^  <  2  x  10  and  A  =  0.0625,  the  core  flow 
structure  Is  non-parallel  and  is  dominated  by  horizontal  intrusions 
flowing  along  each  of  the  two  insulated  horizontal  walls  of  the 
enclosure  and  the  fluid  embraced  by  the  two  horizontal  jets  is 
practically  stagnant  and  thermally  stratified. 
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Sernas  and  Lee  [34,  35]  Investigated  the  heat  transfer  in  a 
rectangular  air  enclosure  both  interferometrically  and  numerically. 
In  [34],  they  investigated  interferometrically  air  in  enclosures 
for  the  range  O.U  A  {  1.0  and  2.64  x  10®  *  Gr^  <  5.45  x  10®. 

Their  enclosures  were  composed  of  two  types  of  horizontal  bounda¬ 
ries.  One  type  was  made  for  isothermal  horizontal  boundary  condi¬ 
tions  and  the  other  type  for  adiabatic  horizontal  boundary 
conditions.  Heat  transfer  characteristics  and  flow  patterns  within 
these  two  types  of  enclosures  were  found  to  be  significantly 
different  for  the  experimental  range  covered.  In  [35],  they  in¬ 
vestigated  the  problem  numerically  for  air  with  0.1  £  A£1.0  and 
10®  £  GrL  £  3  x  107.  When  A  =  0.1,  the  isotherms  and  streamlines 
in  the  isothermal  enclosure  become  very  much  different  from  those 
in  the  adiabatic  enclosure  in  contrast  with  the  case  of  A  *  1.0  and 
0.4  in  which  the  patterns  were  very  similar  to  each  other.  When 
A  =  0.1,  the  flow  pattern  in  the  core  region  is  nearly  parallel  in 
the  adiabatic  enclosure  while  in  the  isothermal  enclosure,  there  is 
almost  no  fluid  motion  in  the  core  except  near  the  end  region. 

Tseng  [36]  presented  some  numerical  results  for  slightly  low 
aspect  ratio  rectangles  in  his  work  on  tilted  rectangular  enclosures 
with  adiabatic,  isothermal  and  perfectly  conductive  connecting  walls 
at  104  £  RaL  £  9  x  107,  0.1  £  Pr  $  100  and  -90°  $  «  $  80°,  For 
RaL  «  10®,  A  =  0.5  and  4  *  0°,  no  appreciable  distinctions  are  ob¬ 
served  in  the  isotherms  and  streamline  distributions  among  those 
three  kinds  of  connecting  walls. 
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Patterson  and  Imberger  [37]  recently  considered  the  problem  of 
transient  natural  convection  in  a  cavity  of  A  £  1.  Scale  analysis 
is  used  to  show  that  a  number  of  initial  flow  types  are  possible, 
collapsing  ultimately  into  two  basic  types  of  steady  flow,  conduc¬ 
tion  dominant  flow  and  convective  flow,  determined  by  the  relative 
value  of  the  non-dimensional  parameters  describing  the  flow.  A 
number  of  numerical  solutions  which  encompass  both  flow  types  are 
obtained,  and  their  relationship  to  the  scale  analysis  is  discussed 
More  recently  Kamotani  et.  al.  [33]  give  a  comprehensive  pic¬ 
ture  of  natural  convection  heat  transfer  in  low  aspect  ratio  enclo¬ 
sures  by  systematically  varying  the  parameters  Ra^,  A  and  Pr.  Among 
the  experimental  results,  secondary  cells  are  present  in  the  end 
region  for  A  =  0.2,  Ra^  =  4.8  x  104  and  Pr  =  965  and  are  shown  to 
influence  the  Nusselt  number  noticeably.  In  such  a  case,  the  core 
flow  pattern  was  not  parallel. 

Most  recently,  Simpkins  and  Dudderar  [38]  experimentally 
studied  convection  in  rectangular  cavity  for  a  variety  of  Rayleigh 
numbers  and  aspect  ratios  for  fluids  with  Prandtl  number  Pr  *  10. 
For  small  aspect  ratios,  e.g.,  for  A  =  0.25,  secondary  rolls  are 
present  and  remain  near  the  outer  edge  of  the  vertical  boundary 
layers  and  the  motion  between  them  becomes  disordered.  There  is 
no  evidence  of  more  than  two  secondary  rolls  appearing  in  the  cavi¬ 
ty.  The  critical  Rayleigh  number  for  secondary  rolls  is  given  by 
(RaH)cA3/4  ^  6.4  x  105  ±  10  %  (1.1) 
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Their  observations  also  suggest  that  the  horizontal  boundary  layers 
have  a  significant  effect  on  the  core  flow  field. 

Shiralkar  and  Tien  [39]  investigated  numerically  heat  transfer 
by  natural  convection  in  a  horizontal  cavity  for  low  aspect  ratios. 
Numerical  results  have  been  obtained  for  Pr  =  1.0,  0.025  -  A  -  1.0 
and  RaH  -  6  x  106  and  for  A  =  0.1,  0.01  iPri  100  and  RaH  <  106. 

For  high  and  moderate  Pr,  the  structure  of  flow  in  the  core  reveals 
essentially  parallel  streamlines.  Nusselt  number  shows  good 
agreement  with  the  prediction  of  Bejan  and  Tien  [30]  at  high  and 
low  RaH-  In  the  intermediate  range,  the  agreement  is  poor.  The 
numerical  results  also  indicate  that  for  low  Pr,  even  at  A  =  0.1, 
at  Ra^  *  10®,  the  core  flow  is  not  fully  developed;  on  the  contrary, 
the  core  contains  developing  velocity  boundary  layers  along  the 
horizontal  walls. 

In  addition,  Jhaveri  and  Rosenberger  [40]  also  studied 
numerically  thermal  convection  in  a  vapor  transport  process  across 
a  horizontal  rectangular  enclosure  for  A  =  0.1,  Sc  *  0.5,  Pr  =  0.7 
a  d  GrH  -  3.3  x  10  ,  which  was  earlier  considered  rather  simply 
by  Klosse  and  Ullersma  [26].  Many  of  the  simplifying  assumptions 
made  in  [26]  were  not  used  In  [40].  The  problem  of  interest  in 
[40]  includes  the  effect  of  mass  flux  as  well  as  heat  flux  and  it 
is  assumed  that  buoyancy  is  due  primarily  to  temperature  differ¬ 
ences.  However,  the  boundary  conditions  differ  from  the  pure 
natural  convection  flows  because  of  the  concentration  difference. 

In  that,  the  horizontal  mass  average  velocity  at  the  Interfaces  is 
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specified  in  terms  of  mass  flux  through  Fick's  law.  The  flow  pro¬ 
files  determined  in  this  way  indicate  the  onset  of  a  recirculation 

4 

cell  at  each  interface  for  Gr^  =  3.3  x  10  .  It  should  be  noted 
that  the  cell  is  due  to  the  coupling  of  mass  flux  at  the  interfaces 
with  the  transport  equations,  otherwise  the  modeled  systems  just 
reduce  to  the  pure  natural  convection  problem. 

1.3  Motivation  of  the  Present  Work 

As  is  pointed  out  in  [1],  in  the  confined  natural  convection 
problems  in  horizontal  cylinders  and  high  aspect  ratio  rectangles, 
two  core  configurations  are  mathematically  possible,  viz.,  rotating 
and  isothermal  or  stagnant  and  stratified.  Analytical  and  numerical 
solutions  first  indicated  the  former  and  it  was  only  after  experi¬ 
ments  indicated  that  the  latter  always  occurred  that  such  results 
were  obtained  theoretically.  Similarly,  in  the  low  aspect  ratio 
rectangular  enclosures,  it  was  shown  experimentally  that  the  flow 
structure  in  the  core  is  not  always  parallel.  Under  certain  condi¬ 
tions,  there  appeared  other  types  of  core  configurations  such  as 
secondary  cells  or  stagnation  flows  in  the  core.  Such  ambiguities 
concerning  the  nature  of  the  flows  are  inherent  in  all  internal  con¬ 
vection  problems  as  has  been  shown  above  and  there  does  not,  as  yet, 
exist  any  way  of  predicting  with  confidence,  a  priori,  what  the  core 
flow  pattern  will  be. 

The  primary  motivation  of  the  present  work  is  to  develop  some 


means  to  predict  the  core  flow  patterns  in  low  aspect  ratio 


rectangular  enclosures. 


1.4  Method  of  Approach 

In  a  low  aspect  ratio  rectangular  enclosure,  the  flow  region 
can  essentially  be  divided  into  two  regions.  One,  the  region  near 
end  walls  which  we  will  call  herein  "end  region"  and  the  other, 

"core  region",  the  region  outside  the  end  region  with  horizontal 
boundaries. 

In  order  to  be  able  to  predict  the  core  flow  pattern,  and  thus, 
the  whole  flow  pattern  correctly,  one  has  to  be  clear  on  what 
physical  mechanisms  govern  the  core  flow.  To  do  this,  one  has  to 
develop  a  mathematical  model  which  properly  represents  the  physics 
of  the  core  flow. 

Since  the  flow  characteristics  may  be  different  and  coupled 
in  each  flow  region,  it  is  not  possible  to  consider  all  the  impor¬ 
tant  physical  mechanisms  in  the  enclosure  by  the  governing  equa¬ 
tions  based  on  either  one  region.  In  addition,  the  equations  that 
describe  the  phenomena  of  interest  herein  are  non-linear  and 
bidirectionally  coupled.  Any  simplifications  must,  therefore,  be 
carefully  made.  In  order  to  ensure  that  an  assumption  made  in  one 
equation  is  consistently  transmitted  to  the  other  equations,  a 
somewhat  formal  procedure  is  employed  herein  based  on  the  method  of 
multiple  scales.  Multiple  scales  are  introduced  to  give  mathemati¬ 
cal  degrees  of  freedom  which  enable  physical  statements  to  be 
properly  made  on  the  important  physical  mechanisms  in  the  enclosure. 
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Analytically  the  method  of  multiple  scales  can  be  used  to  obtain 
single  uniformly  valid  expansions,  in  contrast  with  the  method 
of  matched  asymptotic  expansions  which  yields  two  expansions  that 
must  be  matched  [Nayfeh  (41)]. 

Based  on  the  scaling  analysis  which  gives  degrees  of  freedom, 
proper  force  balances  in  the  basic  dimensionless  equations  can  be 
made  and  the  order  of  magnitude  of  each  physical  terms  can  be 
estimated  with  the  relevant  dimensionless  parameter  as  its  coeffi¬ 
cient.  Then,  by  considering  the  derivatives  with  respect  to  core 
flow  region,  the  equations  which  will  describe  the  core  flow 
characteristics  can  be  extracted  from  the  basic  dimensionless 
equations.  The  core  flow  patterns  are  first  studied  in  a  global 
sense  under  the  implicit  assumption  that  there  exist  no  flow  sub¬ 
regimes,  such  as  secondary  cells.  On  that  basis,  the  geometric 
length  scales  are  the  proper  ones  for  the  core  flow  structure. 

With  the  use  of  the  results  of  the  global  considerations,  detailed 
analysis  is  made  to  determine  the  possibility  of  flow  subregimes. 

In  the  detailed  analysis,  arbitrary  length  scales  replace  the  geo¬ 
metric  length  scales  in  the  dimensionless  equations. 

In  an  attempt  to  predict  the  core  flow  patterns  in  the  detailed 
analysis,  no  criteria  could  be  found  for  the  presence  of  closed 
streamlines  in  the  core  which  would  yield  the  necessary  and  suffi¬ 
cient  conditions  for  the  occurrence  of  secondary  cells.  Instead, 
the  occurrence  of  flow  subregimes  is  considered  to  be  possible 
herein  when  the  arbitrary  horizontal  length  scale,  iQ,  is  less  than 


the  geometric  length  scale,  L,  i.e.,  j~  <  1.  This  statement  may 
not  be  the  sufficient  condition,  but,  at  least  it  can  be  considered 
to  be  a  necessary  condition  for  the  occurrence  of  flow  subregimes. 
The  results  obtained  on  this  premise  are  compared  with  existing 
experimental  data. 


CHAPTER  II 


MATHEMATICAL  FORMULATION  OF  THE  PROBLEM 


The  problem  to  be  studied  herein  is  the  convective  motion  of  a 
fluid  driven  by  the  imposed  temperature  difference  between  the  end 
walls  within  a  two-dimensional  rectangular  enclosure  of  length  L 
and  height  H.  The  schematic  diagram  of  the  system  is  shown  in  Fig. 

1. 

The  flow  field  will  be  considered  to  be  steady,  two-dimensional, 
laminar  and  quasi-incompressible.  By  quasi-incompressibility  it  is 
meant  that  the  variations  in  density  due  to  temperature  differences 
are  neglected  except  where  they  modify  the  body  force  terms  in  the 
momentum  equations.  This  is  generally  referred  to  as  Boussinesq 
approximation  and  has  been  discussed  in  a  formal  way  by  Ostrach  [42} 

2.1  Normalization  of  the  Basic  Equations 


The  basic  equations  of  vorticity  and  energy  transport  for  two- 
dimensional  flow  can  be  written  as 


n 
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(2.1) 
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Figure  1.  Schematic  Diagram  of  the  System 
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where  5  denotes  vorticity,  v  the  stream  function  and  T  temperature. 
Also  g  is  the  acceleration  due  to  gravity  and  v,  6  and  o,  the 
kinematic  viscosity,  coefficient  of  thermal  expansion  and  thermal 
diffusivity. 

The  viscous  dissipation  was  neglected  in  the  energy  equation 
(2.3)  due  mainly  to  the  very  low  velocities  present  in  natural 
convection  flows  of  the  type  considered  herein.  The  variations  of 
the  dynamic  viscosity  v  and  thermal  conductivity  k  were  also 
neglected  as  a  result  of  the  quasi-incompressiblity  of  the  fluid. 

The  corresponding  boundary  conditions  are 

£  =  s  0  -  T  =  Tw<  Te  *  X  -  o,  L  (2.4a) 

it  *  -  o  j  *  O  at  Y  *  H  (2.4b) 


The  basic  equations  (2.1)  '>■  (2.3)  are  made  dimensionless  by 
the  following  definitions. 


x 


X 

L 
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-T 


(2.5) 


In  (2.5),  the  characteristic  stream  function  vR  is  to  be  specified 

later  in  the  analysis  in  accordance  with  the  physics  of  the  problem 

2 

The  characteristic  vorticity  =R  is  represented  as  vR/t  ,  where  i  is 
a  length  scale  that  will  also  be  specified  in  the  course  of  analy¬ 
sis.  The  purpose  of  introducing  the  unspecified  length  scale  ; 
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(instead  of  using  the  specific  geometric  length  scale)  is  to 
normalize  the  vorticity  properly  according  to  the  physical  situation 
of  the  system. 

Substitution  of  the  definitions  (2.5)  into  eqs.  (2.1)  n,  (2.3) 
yields 


38$-  + 

£i-  ) 

3?*  / 

(2.6) 

5  =  (  A*  w*  +  ^  ) 

(2.7) 

2LZ11  -  «L  /  S- +  3flL\ 

'dcx.,))  v  "  -it1  T  ) 

(2.8) 

u 

where  aT  =  TH-TC  and  A  represents  the  aspect  ratio,  p 

The  corresponding  boundary  conditions  (2.4)  become 

f  -  T*  *  9  -  0 

at  X*  0 ,  1 

(2.9a) 

-<5r 

It 
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II 

'C. 

0,  1 

(2.9b) 

One  property  to  be  pointed  out  herein  is  that  for  adiabatic 
(or  insulated)  horizontal  boundaries,  the  equations  and  boundary 
conditions  have  a  symmetry  property  involving  reflections  about 
x  *  1/2,  y  =  1/2  of  the  cavity.  This  was  termed  the  centro- symmetry 
property  of  the  equations  and  boundary  conditions  by  Gill  [7],  and 
can  be  represented  as 


L 
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f<*.  V  -  ftl-X.  1-)) 

(2.10a) 

91*-.})  -  1-  i-f ) 

(2.10b) 

2.2  Working  Form  of  the  Equations 

2 

In  the  dimensionless  equations  (2.6)  (2.8),  when  A  <<  1,  the 

x-direction  diffusion  terms  can  be  neglected  compared  to  the  y- 
direction  diffusion  terms.  Since  the  x-direction  diffusion  terms 
are  the  highest  order  x-derivatives,  if  they  are  neglected,  the  end 
wall  boundary  conditions  can  not  be  satisfied.  Therefore,  the 
equations  without  the  horizontal  diffusion  terms  in  Eqs.  (2.6)  -v 
(2.8)  are  only  valid  away  from  the  ends,  i.e.,  in  the  core.  On  the 
other  hand,  the  equations  which  represent  the  flow  characteristics 
in  the  end  region  could  be  obtained  by  properly  stretching  the  x- 
coordinate  in  the  Eqs.  (2.6)  M2. 8). 

In  order  to  obtain  insights  into  the  nature  of  flow  character¬ 
istics,  the  unspecified  characteristic  quantities  such  as  Vp  and 
*  in  the  Eqs.  (2.6)  ^  (2.8)  first  should  be  determined  to  see  what 
the  relevant  dimensionless  parameters  are.  These  quantities  may, 
in  principle,  be  determined  by  considering  the  flow  mechanisms  in 
either  region.  However,  as  the  two  regions  are  closely  coupled  and 
this  coupling  affects  much  of  the  flow  structure,  it  is  necessary 


that  in  determining  the  characteristic  quantities,  both  flow  regions 
and  their  coupling  be  considered  in  some  way.  Otherwise,  some 
important  physical  mechanisms  may  be  overlooked  by  neglecting  the 
coupling  effects  on  the  flow  characteristics.  A  mathematically 
formal  procedure  is,  therefore,  employed  based  on  the  method  of 
multiple  scales  through  which  both  the  end  and  core  flow  regions 
together  with  the  interacting  region  are  explicitly  identified  in 
the  basic  dimensionless  equations. 

Here  the  physical  domain  of  interest  is  confined  and  finite 
and  the  singular  behaviour  of  the  dimensionless  equations  are 
expected  near  the  end  walls.  Thus  the  multiple  scales  are  intro¬ 
duced  in  the  following  way. 

3  •  *  ,  7  -  -t  <2-u> 

and  the  derivatives  are 

=  J—  X-  +  3- 

'bX,  £*  07  T  'iy  (2.12) 

V  _  I  .  X  .  V 

t?  7?  4?  +  X  w  + 

where  the  derivatives  with  respect  to  ;  and  n  represent  the  core  and 
end  region  characteristics  respectively  and  ex  is  the  small,  stretch¬ 
ing  parameter. 

Introducing  the  derivatives  of  multiple  scales  (2.12)  into  the 
dimensionless  equations  (2.6)  n.  (2.6),  we  then  obtain 
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(2.13) 

(2.14) 


7*0  (2.15) 


Those  are  the  working  form  of  the  basic  dimensionless  equations. 

In  the  above  equations,  it  is  possible  that  the  extent  of  the 
end  region  may  be  different  for  the  temperature  and  vorticity  dis¬ 
tributions.  This  introduces  some  ambiguity  in  cx  used  in  the 
coordinate  stretching.  It  thus  should  be  kept  in  mind  that  n  is 
not  a  normalized  variable  that  always  yields  the  proper  normaliza¬ 
tion  of  each  term  in  the  end  region  and  must,  therefore,  be  care¬ 
fully  treated.  Some  modification  will  be  made  if  necessary  for  the 
proper  normalization  in  accordance  with  the  physical  situation 
concerned  therein. 

In  the  following  analysis,  the  dimensionless  groups  are  first 
to  be  determined  according  to  three  different  cases  of  force 
balances  in  Eqs.  (2.13)  ^  (2.15).  Then  the  core  flow  characteris¬ 
tics  are  going  to  be  studied  from  a  global  view  by  first  noting  the 
character  of  the  dimensionless  groups,  then  observing  their  effects 
on  the  equations  as  they  are  allowed  to  take  on  limiting  values. 
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Finally,  the  core  flow  pattern  will  be  examined  in  detail  with  the 
information  obtained  from  the  global  analysis. 
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CHAPTER  III 


GLOBAL  VIEW  OF  THE  FLOW  CHARACTERISTICS 


CASE  I :  Balance  Between  Buoyancy  and  Viscous  Forces  in  the  Core 

For  a  given  fluid  and  geometry,  when  there  exist  no  distinct 
thermal  layers  at  the  end  walls  within  which  most  of  the  tempera¬ 
ture  change  occurs,  the  imposed  temperature  gradient  will  be  felt 
across  the  entire  configuration.  Since  the  driving  force  of  the 
fluid  motion  concerned  herein  is  the  buoyancy  force  induced  only  by 
the  temperature  gradient,  in  this  situation  for  low  aspect  ratio 
enclosures,  i.e.,  for  A  «  1,  the  driving  force  will  be  the  buoy¬ 
ancy  force  in  the  core.  We  therefore  study  this  situation  first 
considering  the  balance  between  buoyancy  and  viscous  forces  in  the 
core.  The  explicit  conditions  for  each  situation  concerned  herein 
will  be  delineated  in  the  course  of  analysis. 

3.1.1  Basic  Equations 


From  Eq.  (2.13),  the  balance  between  buoyancy  and  viscous 
force  in  the  core  can  be  represented  as 


ggAT*  H  vL_ 

V 


(3.1) 


Considering  that  the  buoyancy  force  in  the  core  drives  the 
core  flow  and  that  there  is  not  supposed  to  exist  any  particular 
region  in  the  core  within  which  the  viscous  forces  are  so  concen- 
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centrated  as  to  balance  the  core  flow,  it  is  appropriate  to 
represent  the  characteristic  vorticity  ,  by  specifying  the 
characteristic  length  i  as  H,  as 

t  -v  H  (3.2a) 


and 


-R 


-v  R 
hT 


From  (3.1)  and  (3.2),  we  obtain 

7  -v  89-TH4 
*R  vL 


(3.2b) 


(3.3) 


One  more  balance  is  required  here  to  determine  the  stretching 
parameter,  e ,  from  which  the  end  region  flow  equations  and  its 
characteristic  length  can  be  determined.  However,  at  this  moment, 
that  is  not  clear  before  we  know  the  physics  of  the  core  flow 
properly.  Thus  this  will  be  considered  later  after  studying  the 
core  flow  characteristics. 

Substituting  (3.3)  into  Eqs.  (2.13)  ^  (2.15)  and  considering 
the  terms  with  z,  and  y  derivatives,  the  equations  which  will 
describe  the  core  flow  characteristics  can  be  written  as 


5 


(3.4) 

(3.5) 

(3.6) 
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For  the  situation  described  above,  it  is  now  evident  that  our 

2 

analysis  is  valid  under  the  condition  Gr^A  £  1.  Beyond  this 

2 

condition,  i.e.,  for  Gr^A  >>  1,  inertia  becomes  dominant  in  the 
core,  and  thus,  different  scalings  are  required  based  on  differ¬ 
ent  force  balances  for  the  analysis  of  those  situations. 

Based  on  the  dimensionless  parameters  of  Eqs.  (3.4)  (3.6), 

we  are  going  to  study  the  core  flow  characteristics.  After  that, 
some  consideration  will  be  given  on  the  end  region  based  on  the 
core  configurations. 

3.1.2  Core  Flow  Characteristics 

[1]  GrHA2  «  1,  A2  «  1 

Under  this  condition,  the  core  flow  equations  (3.4)  ^  (3.6) 
reduce  to 


o  =  +  XL 

u  *5  'iT 
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(3.7) 

(3.8) 

(3.9) 


Because  Pr  appears  as  a  parameter  in  Eq.  (3.9),  its  effect  must  be 
analyzed.  We  thus  examine  the  flow  characteristics  according  to  Pr. 
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(i)  Pr  £  1  -  RaHA2  «  1 

2 

When  Pr  x  1  such  that  Ra^A  <<  1,  from  Eq.  (3.9) 

0  -  (3.10) 

Considering  the  adiabatic  horizontal  boundary  conditions  (2. 9b) 

6  =  0(5)  (3-H) 

%  fa 

Since  the  buoyancy  term  in  the  core,  in  Eq.  (3.7)  depends 

on  the  variation  of  c,  its  dependence  on  c  should  be  known  to  obtain 
the  core  streamline  profiles.  From  (3.11),  considering  that  the 
heat  flow  is  one-dimensional  in  a  purely  conductive  regime, 
i.e.,  one-dimensional  conduction  heat  transfer,  9  should  be  linear 
with  respect  to  c  as 

6  =  K,  J  +  K*  (3.12) 

where  k^  and  k£  are  arbitrary  constants. 

From  Eqs.  (3.7),  (3.8)  and  (3.12),  we  then  have 

Considering  the  boundary  conditions  (2.9b),  from  (3.13)  can 
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be  represented 


r  *&»*(»-/>* 


(3.14) 


From  (3.12)  and  (3.14),  it  is  seen  that  the  temperature 
gradient  in  the  core  which  provides  the  driving  force  for  the  core 
flow  is  linear  and  the  corresponding  core  flow  pattern  is  parallel. 

(ii )  Pr  »  1  -  RaHA2  1 

2 

When  Pr  is  large  such  that  Ra^A  a.  1,  Eq.  (3.9)  reduces  to 


2UL12  =  ££- 

*(5,  V 


(3.15) 


Since  the  flow  characteristics  from  Eq.  (3.15)  are  not  readily 
perceivable,  we  try  to  get  information  from  an  integral  form  of  Eq. 
(3.15)  with  the  boundary  conditions  at  the  horizontal  boundaries. 

If  we  rewrite  Eq.  (3.15) 


(3.16) 


Integration  of  Eq.  (3.16)  with  respect  to  y  from  y  =  0  to  1 
yields 


(3.17) 


i 


e1 


II 
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By  changing  the  order  of  integration  and  differentiation  and 
applying  the  Leibnitz's  rule  in  the  first,  two  expressions  respect¬ 
ively  on  the  left  side,  Eq.  (3.17)  can  be  rewritten  as 


The  left  side  is  ingetrated  by  parts  and  we  then  have 


(3.19) 


Applying  the  boundary  conditions  (2.9b)  for  v  ,  Eq.  (3.19)  becomes 


(3.20) 


For  adiabatic  horizontal  boundaries  or  for  the  equal  heat  flux  at 
both  upper  and  lower  horizontal  boundaries,  i.e.,  for 


or 


(3.21) 

(3.22) 

(3.23) 


The  situations  for  which  either  Eq.  (3.22)  or  Eq.  (3.23)  can 
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be  satisfied  are  examined  in  Appendix  A.  From  the  considerations 
made  therein,  it  is  shown  that  the  only  possible  configuration 
is  that  the  temperature  distribution  in  the  core  is  linear  and 
stratified,  while  the  core  flow  pattern  is  still  parallel,  i.e., 
e  and  v  are  represented  as 

B  *  K3  5  +  9lV  (3.24) 

and  f  =  (3.25) 

where  is  an  arbitrary  constant  and  g(y)  is  an  arbitrary  function 
of  y. 

Substitution  of  (3.24)  and  (3.25)  into  Eqs.  (3.7)  and  (3.8) 
yields 

h "'('})  =  /  (3.26) 

and  from  the  boundary  condition  (2.9b),  h(y)  can  then  be  represen¬ 
ted  as 

=  (3.27) 

[2]  GrHA2  1,  A2  «  1 

Under  this  condition,  the  core  flow  equations  (3.4)  n. 
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(3.6)  reduce  to 


W-V)  _  39  . 

-as  T 

u_ 

(3.28) 

i 

ii 

>*> 

(3.29) 

P-M‘  liLIL  - 

3(5,  J) 

re 

(3.30) 

As  in  the  preceding  analysis,  the  flow  characteristics  are  being 
examined  according  to  Pr. 

(i)  Pr  1  -  RaHA2  <<-  1 

o 

When  Pr  is  very  small  such  that  Ra^A  <<  1,  from  Eq.  (3.9) 

0  -  (3.31) 

This  equation  is  identical  to  Eq.  (3.10)  which  was  previously 
treated.  Referring  the  result  in  (3.12),  e  can  be  represented  as 

B  =  M  +  K  (3-32) 


By  the  same  reasoning  as  in  the  case  (ii)  of  the  preceding 
analysis,  if  we  rewrite  and  integrate  Eq.  (3.28)  over  y  from  0  to 
1  we  have 


(3.33) 


30 


Integrating  the  first  expression  on  the  left  side  by  parts, 


(3.34) 


Substituting  Eq.  (3.8)  into  Eq.  (3.34)  and  applying  the  boundary 
conditions  (2.9b)  for  ,  we  then  have 


jL. 

^5 


£L 

-if1 


A*  + 


(3.35) 


Since  the  flow  characteristics  contained  in  Eq.  (3.35)  are  not 
readily  perceivable,  we  attempt  to  obtain  information  by  examining 
the  dependence  of  c  on  the  arguements  ?  and  y  in  Eq.  (3.35). 

Considering  the  boundary  conditions  (2.9b),  c  can  be  represen¬ 
ted  as 


ff5,n  =  (3.36) 

From  the  symmetry  property  of  ^  in  (2.10a),  we  have 


kc  5) 

-  hci'1) 

(3.37a) 

and 

=  *</- f ) 

(3.37b) 

Substitution  of  (3.32)  and  (3.36)  into  Eq.  (3.35)  yields 


f 
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+  hli)[  f '>'!){}  -  h  (3.38) 

Carrying  out  the  differentiation  in  the  first  expression  on  the  left 
side,  we  find 


i 


L 


2hci)k'(f)f&lV)f'(t)(jz  +  hot)  [}'"(*)!']  *  K,  (3.39) 

In  Eq.  (3.39),  considering  the  boundary  conditions  and  the 
symmetry  property  of  y  in  (2.9b)  and  (3.37),  the  function  g(y) 
is  an  even  function  about  the  point  y  =  $.  The  derivative  g'"(y) 
is  an  odd  function  about  y  =  J,  and  thus  the  product,  g(y)g"'(y), 
becomes  an  odd  function  about  y  =  i.  Then  the  integral 
;og(y)g"'  (y)  dy  becomes  zero  and  Eq.  (3.39)  reduced  to 


k(V  (?'"< -  «, 

(3.40) 

Thus 

hit)  =  *jp°—  -  con$t.  - 

(3.41) 

Then  from  (3.36)  and  (3.41),  we  have 

f 

(3.42) 

Substituting  (3.32)  and  (3.42)  into  Eqs.  (3.28)  and  (3.29),  from 
the  boundary  conditions  (2.9b)  g(y)  can  be  represented  as 


* 


(3.43) 
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It  is  thus  seen  from  (3.32)  and  (3.42)  that  the  core  flow 
pattern  is  parallel  and  the  core  temperature  distribution  is  linear. 

(ii)  Pr  ^  1  Ra^A2  -v  l 

2 

When  Pr  is  of  order  1  such  that  Ra^A  ^  1 ,  Eq.  (3.9)  can 
be  written  as 


=  XI 

mi,})  ~ 


(3.44) 


This  equation  is  identical  to  Eq.  (3.15)  which  was  previously 
treated.  Referring  the  results  indicated  in  Appendix  B,  ^  and  e 
can  be  represented  as 

f  =  (3.45) 

B  -  fCS)  t  ?l})  (3.46) 

Substitution  of  (3.45)  and  (3.46)  into  Eqs.  (3.28)  and  (3.29)  yields 

fit)  =  fVC»)  =  K¥  (3.47) 

where  is  an  arbitrary  constant. 

From  (3.46)  and  (3.47),  6  is  represented  as 

9  *  +■  (3.48) 

From  (3.47)  and  the  boundary  conditions  (2.9b),  $  then  becomes 

ts  fo-' f 


(3.49) 


r 
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Thus  from  (3.48)  and  (3.49),  it  is  seen  that  the  core  flow  pattern 
is  still  parallel  in  this  case  while  the  core  temperature  distri¬ 
bution  is  linear  and  stratified.  These  are  the  same  results  as 

2 

those  in  the  case  (ii)  of  the  preceding  analysis  for  Gr^A  «  1, 

2 

Ra^A  -v  1  . 

2 

(iii )  Pr  »  1  -*•  Ra^A  »  1 


In  this  situation,  from  Eq.  (3.6),  heat  transfer  by  conduc¬ 
tion  is  negligible  and  convection  is  dominant  in  the  core.  Conduc¬ 
tion  of  heat  from  the  hot  wall  will  thus  be  transferred  by  convection 
across  the  cavity  to  the  opposite  end.  This  implies  the  existence 
of  thermal  boundary  layers  in  the  end  region  within  which  most  of 
the  temperature  drop  occurs  and  there  no  longer  exists  driving 
force  in  the  core.  Instead  the  driving  force  will  come  out  from 
the  end  region.  For  the  analysis  of  this  flow  regime,  different 
scalings  based  on  the  proper  balances  in  the  end  region  are  needed. 
This  situation  will  be  considered  as  Case  III. 

[3]  GrHA2  »  1,  A2  «  1 


In  this  situation,  from  Eq.  (3.4)  inertia  becomes  dominant  in 
the  core.  Thus  the  balance  made  between  buoyancy  and  viscous  forces 
in  the  core  becomes  inappropriate  to  this  situation.  The  flow 
inertia  can  be  produced  either  by  the  driving  force  in  the  end 
region  or  by  the  driving  force  in  the  core.  The  former  case  will 


be  considered  as  Case  III  and  the  latter  case  will  be  considered  as 


Case  II.  The  explicit  conditions  for  each  situation  will  be 
delineated  in  the  course  of  analysis. 

3.1.3  Determination  of  Stretching  Parameter  ex 


In  the  preceding  analysis,  it  was  shown  that  the  temperature 
distribution  in  the  core  is  either  linear  or  linear  and  stratified 
and  its  flow  pattern  in  the  core  is  parallel.  In  such  a  flow 
regime,  the  driving  force  for  the  core  flow  results  from  the  linear 
temperature  gradients  in  the  core  and  the  end  region  is  supposed 
to  have  little  effect  on  the  core  flow  except  simply  to  turn  the 
core  flow  to  conserve  the  mass.  Thus  the  proper  force  balance  in 
the  end  region  that  is  needed  to  determine  ex  would  be  that  between 
the  viscous  diffusion  along  the  end  walls  and  that  along  the 
horizontal  boundaries.  From  Eq.  (2.13),  this  balance  can  be 
represented  as 


and  thus,  ex  <v  A  (3.51) 

From  Eq.  (2.11),  the  variable  x  in  the  end  region  is  of  order 

«x 

0(— )  and  it  must  be  stretched  to  make  n  of  unit  order  of  magni¬ 

tude  by  means  of  the  stretching  parameter,  e  ,  so  that 


Then,  from  (3.51)  and  (3.52),  the  end  region  characteristic  length 


6x  becomes  of  order 


6  ■v  e  •  L  'v.  H 

x  x 


(3.53) 


Substitution  of  (3.51)  and  (3.3)  into  Eqs.  (2.13)  <v  (3.15) 
yields 


Cr  fi2LU±  +  <*1*21112  =  afi_  JL + 


+  **  +  A 


^  ^  +  2A  vfai  +  +  ir  ) 


>5* 


(3.54) 

(3.55) 


p,r  fUSxtl+Pti.  /t2ML  =  A3&-tiAL3S-  +/^fi-  +  A-2^r  (3  56) 

£ ^ A  tnTfj *  */ < hT  *  s  v  r  *  -ar  ar  50 ' 


57u)  ' n  *y 


>r 


Considering  the  terms  with  and  y  derivatives  in  the  Eqs. 
(3.54)  -v  (3.56),  the  equations  which  will  describe  the  end  flow 
characteristics  can  be  extracted.  Before  doing  that,  one  thing 
that  should  be  noted  is  that  in  Eq.  (3.54),  by  a  coordinate  stretch 


ing  with  a  small  parameter  A,  the  buoyancy  term  in  the  end  region, 


1  80  9  ^ 

j  -j^,  becomes  dominant  relative  to  the  core  buoyance  term, 


which  provides  the  driving  force  to  the  core  flow.  Physically  this 
is  unreasonable.  This  is  actually  due  to  the  fact  that  the  tem¬ 
perature  difference  in  the  end  region,  which  is  a  small  part  of 
the  total  temperature  change  aT,  is  not  properly  normalized  by  a 
simple  stretching  of  the  independent  variable  x.  This  is  one  of 
the  difficulties  encountered  due  to  the  coupling  of  the  governing 
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equations,  because  in  the  coupled  equations  such  as  Eqs.  (2.13)  and 
(2.15)  herein,  the  two  dependent  variables  e  and  •;  do  not  behave 
similarly  in  the  region  of  non-uniformity,  and  thus  the  stretching 
parameter,  ex,  which  is  determined  from  the  consideration  of  either 
one  dependent  variable,  does  not  normalize  the  terms  with  the  other 
variable  properly  in  the  region  of  non-uniformity.  Some  modifica¬ 
tion,  therefore,  should  be  made  to  yield  the  proper  normalization  of 
the  relevant  equations  if  necessary. 

In  Eq.  (3.54),  if  we  separate  the  temperature  distribution  into 
two  parts  as 


%.  V  =  *  9CU.  })  (3-57) 


where  the  subscripts  c  and  e  represent  the  core  and  end  region, 
respectively,  since 


and 


we  find 


6  - 


ST-Til 

AT 


£L 

AT 


V 


I 


A 


*  IT.-JDc 


A  AT 
AT 


(3.58) 

(3.59) 


(3.60) 


Thus,  for  the  proper  normalization  of  the  end  region  buoyancy  term, 
we  replace  ~  by  A  in  Eq.  (3.54). 
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Then,  considering  the  terms  with  n  and  y  derivatives,  the  end 
region  flow  equations  can  be  written  as 


2LU2  -  U.  +  :£i_  +  2lJ_ 

5  =  _  (  iL  +  iL\ 

?  '  vy  r  'iy  ) 


R-Cfr  a  21L1L  =  t  A: 


(3.61) 

(3.62) 

(3.63) 


The  method  of  multiple  scales  has  yielded  Eqs.  (3.54)  -v  (3.56) 
which  indicate  explicitly  what  equations  are  applicable  in  each  re¬ 
gion  and  in  the  interaction  region.  From  those  equations,  two  sets 
of  equations,  one  for  the  core  flow  and  the  other  for  the  end  region 
flow,  are  extracted.  Eqs.  (3.4)  -v  (3.6)  for  the  core  flow  and  Eqs. 
(3.61)  -v  (3.63)  for  the  end  region  flow  are  identical  to  those 
obtained  by  Cormack  et.  al.  [27]  for  A  -*  0,  Gr^  fixed  while  present 
equations  are  valid  for  A  <<  1,  GrHA  £  1  and  RaHA  *  1.  However, 
the  buoyancy  term  was  not  properly  normalized  in  [27]  so  that  the 
coefficient  in  their  equations  differ  from  that  in  Eq.  (3.61).  It 
should  be  noted  that  the  buoyancy  does  contribute  to  the  flow  in  the 
end  region  in  contradistinction  to  the  description  given  by  Cormack 
et.  al. 

Since  our  primary  concern  herein  is  general  flow  patterns  based 
on  the  core  configurations,  no  attempt  will  be  made  to  obtain  the 
solution  of  the  basic  equations  at  this  time. 
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3.2.1  Basic  Equations 

In  Case  I,  the  analysis  was  made  based  on  the  balance  between 

buoyancy  and  viscous  forces  in  the  core.  The  balance  made  therein 

2 

was  shown  to  be  inappropriate  when  Gr^A  >>  1,  because  the  inertia 
becomes  dominant  in  the  core  under  that  condition.  We  now  study 
the  flow  characteristics  considering  the  balance  between  buoyancy 
and  inertia  forces  in  the  core. 

From  Eq.  (2.13),  the  balance  between  buoyancy  and  inertia 
forces  in  the  core  can  be  represented  as 


(3.64) 


Specifying  t  as  H  by  the  same  reasoning  in  (3.2),  vR  becomes 


vR  ^  (sgiTH3)* 


(3.65) 


Here  we  also  leave  the  determination  of  the  stretching  parameter, 
cx>  to  later,  after  studying  the  core  flow  characteristics. 

Substituting  (3.65)  into  Eqs.  (2.13)  n.  (2.15)  and  considering 
the  terms  with  i  and  y  derivatives,  the  core  flow  equations  can  be 
written  as 


Mil  =  m.  + 

W,'})  TS 


( yf  TlL.  +  j£L.  ) 
^  'bV  ' 


(3.66) 
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(3.67) 


mJi  =  /  —  f  +  2&- ) 

iPrR^^  K  w  J 


(3.68) 


For  the  situation  Gr^A  >>  1  concerned  herein,  Pr  appears  as 
a  parameter  in  Eq.  (3.68).  Because  consideration  is  being  given 

to  the  situation  in  which  the  driving  force  exists  in  the  core,  the 

2 

analysis  will  be  valid  within  the  parameteric  range  of  PrRaHA  f  1 

2  1 

which  is  possible  only  for  Pr  ^ 


GrHA 


,  i .e. ,  for  low  Pr.  But  it 


will  be  shown  later  in  the  analysis  that  the  validity  of  analysis 

2  2  1 

is  confined  to  the  parametric  range  of  Gr  A  •>  1,  fr  ••  - * 

2  GrHA 

so  that  PrRa^  <<  1.  In  addition,  one  thing  to  be  noted 

2 

is  that  for  GrHA  >>  1,  viscous  diffusion  terms  would  seem  to  be 

negligible  compared  to  the  other  terms  in  Eq.  (3.66).  However, 

since  these  terms  are  the  derivatives  of  the  highest  order  in  the 

equations,  Eq.  (3.66)  is  singular  and  this  implies  that  for  negli- 

2  2 

gible  horizontal  diffusion  term,  A  — *  ,  in  the  core  for  A  <<  1, 

3  r 

there  exists  a  thin  layer  very  near  the  horizontal  boundaries 

32£ 

within  which  the  vertical  viscous  diffusion  term,  — |  ,  is  impor- 

sy 

tant.  The  effect  of  viscous  diffusion  is  not  regarded  to  be  so 
large  as  to  affect  the  core  flow  structure  in  this  inertia  impor¬ 
tant  flow  regime  concerned  herein. 
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3.2.2  Core  Flow  Characteristics 

For  Sr^A^  >>  1,  Eq.  (3.66)  reduces  to 


liLli  =  2L 
Li,})  'i  5 


(3.69) 


2  1  2 

For  Pr  £  - *•  ,  PrRauA  acts  as  a  parameter  in  Eq.  (3.68).  We 

Or  A2  H 

thus  examine  the  core  configurations  according  to  this  parameter, 
(i )  PrRaHA2  <<  1 ,  ' '  <<  1 


2  *  2 
When  Pr  <<  — ^  ,  i.e.,  Pr  is  very  small  so  that  PrRaHA 

GrHA 

<<  1,  Eq.  (3.18)  reduces  to 


0  = 


XL 


(3.70) 


This  equation  is  identical  to  Eqs.  (3.10)  and  (3.31)  which  were 
previously  treated.  Referring  the  result  in  (3.18),  e  can  be 
represented  as 


&  =  K,  3  +  k,  (3.71) 

Then,  from  Eq.  (3.19)  and  (3.71) 

ULLL  -  u  -  const.  (3.72) 

In  order  that  Eq.  (3.72)  is  valid,  i|<  should  be  a  function  of  both 
variables  as 


f  =■  Yu.  7) 


(3.73) 
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This  implies  that  the  core  flow  pattern  becomes  non-parallel  when 
the  inertia  is  important  in  the  core,  although  the  flow  is  driven 
by  the  buoyancy  force  in  the  core. 

(ii )  PrRaHA2  n,  1,  A2  «  1 


In  this  case,  Eq.  (3.68)  reduces  to 


2UL11  -  VL 

at},*;  " 


(3.74) 


This  equation  is  identical  to  Eq.  (3.15)  in  Case  I  which  was 
previously  treated.  Following  the  same  procedure  in  Case  I,  for 
adiabatic  horizontal  boundaries,  from  Eqs.  (3.22)  and  (3.23)  we  have 


or 


o 

}'( +  • 
0 


The  situations  for  which  either  Eq.  (3.75)  or  (3.76)  can  be 

satisfied  are  examined  in  Appendix  B.  From  the  considerations 

made  therein,  it  was  found  that  the  driving  force  for  this  case 

could  not  occur  in  the  core.  Thus,  in  this  situation  the  balance 

based  on  the  driving  force  in  the  core  is  inappropriate  and  the 

2 

analysis  based  on  that  balance  will  be  valid  only  for  PrRaHA  <<  1. 
2  1 

Actually  for  Pr  « - (i.e.,  for  Pr  <<  1),  the  condition 

2  Gr^A  2 

PrRaHA  l  is  possible  only  when  RaHA  »  1  in  which  case  the 

driving  force  was  mentioned  to  exist  in  the  end  region  in  Case  I. 
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This  situation  will,  therefore,  be  considered  as  Case  III. 

(i i i )  PrRaHA2  >>  1 

In  this  case,  convection  becomes  aominant  in  the  core  in  Eq. 

(3.68),  thus  there  no  longer  exists  the  driving  force  in  the  core. 

2 

This  is  also  the  case  of  RauA  >>  1  and  will  be  considered  as  Case 

in. 


3.2.3  Determination  of 

In  the  above,  although  the  core  flow  is  driven  by  the  buoyancy 
force  in  the  core,  the  core  flow  pattern  was  shown  to  be  non-paral¬ 
lel.  In  this  situation,  the  proper  force  balance  for  the  end 
region,  from  which  the  stretching  parameter,  c  ,  is  also  Determined, 
is  not  readily  apparent.  However,  considering  the  relative  domin¬ 
ance  of  inertia  with  respect  to  viscous  friction  in  the  core  and 
the  consequent  non-parallel  core  flow  structure,  two  possible 
mechanisms  can  be  inferred:  The  first  one  is  that  the  inertia  in 
the  end  region  can  be  as  influential  on  the  core  flow  structure  as 
the  inertia  in  the  core.  The  second  one  is  that  the  inertia  in  the 
core  can  be  damped  by  the  viscous  diffusion  in  the  end.  Based  on 
these  two  inferrences,  we  thus  make  the  corresponding  force  bal¬ 
ances  for  the  end  region  and  to  see  what  the  results  would  imply. 


For  the  first  situation,  the  balance  between  the  inertia  in 
the  end  and  the  inertia  in  the  core  is  appropriate.  From  Eq.  (2.13), 


43 


this  balance  can  be  represented  as 

~  %  1  (3.77) 

x 

Thus,  cx  "v  1  (3.78) 

From  (3.52)  and  (3.78),  6^  becomes  of  order 

ex  •  L  L  (3.79) 


This  tells  us  that  there  is  no  difference  in  horizontal  length 
scale  between  the  core  and  end  region,  which  means  that  no  clear 
distinction  could  be  made  between  the  two  regions.  In  addition, 
since  ex  1,  there  is  no  distinction  between  the  equations  with  the 
variables  (r,,  y)  and  (c,  y),  i.e.,  in  the  end  ar,d  core  regions. 

The  core  flow  picture  in  this  situation  will  thus  be  valid  through¬ 
out  the  whole  cavity  except  the  very  thin  layer  adjacent  to  the 

solid  boundaries  which  is  supposed  to  exist  due  to  the  singular  be- 

? 

haviour  of  the  viscous  diffusion  terms  for  Gr^A  >>  1  in  Eq.  (3.66). 
Consequently  the  vertical  flow  (or  turning  flow)  in  the  cavity  is 
not  necessarily  confined  to  the  region  near  the  end  walls  in  this 
case  as  in  the  case  of  parallel  core  flow.  Instead  it  may  spread 
over  the  entire  configuration. 

For  the  second  situation,  the  balance  between  the  viscous 
diffusion  in  the  end  and  the  inertia  in  the  core  is  appropriate. 

From  Eq.  (2.13),  this  balance  can  be  represented  as 


uL 


'V  1 


(3.80) 
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From  (3.65)  and  (3.80),  becomes 
X  (GrHA2)* 

From  (3.52),  & then  becomes  of  order 


(3.81) 


(3.82) 


J<  ^  .  (3.83) 

H  1 

Substituting  (3.65)  and  (3.81)  into  Eqs.  (2.13)  ^  (2.15)  and 
considering  the  terms  with  n  and  y  derivatives,  the  end  flow  equa¬ 
tions  can  be  extracted.  However,  in  this  situation  the  end  region 

1  3  P 

buoyancy  term,  — —  ,  which  is  of  order 

ex  an 

i_  3^  v  (GrHft2)i  16  (3.84) 

ex  an  ^  A  3n 

becomes  dominant  relative  to  the  core  buoyancy  term,  . 

As  in  the  analysis  of  Case  I,  even  though  we  renormalize  the 
end  region  buoyancy  term  with  AaT,  the  proper  end  region  characteris¬ 
tic  temperature  difference,  the  modified  end  region  buoyancy  term, 

(Gr  A^)*  —  *  becomes  greater  than  that  in  the  core  for  the 

condition  Gr^  >>  1  herein.  This  implies  that  temperature  drops 
greater  than  AaT  occurs  in  the  end  region  ,  and  that  will  be  possible 
when  rather  steeper  non-linear  temperature  gradient  exists  in  the 
end.  In  this  case,  considering  the  fact  that  -g  <  1  in  (3.83),  the 
core  flow  will  be  rather  driven  by  the  buoyancy  in  the  end  region 


i 
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than  by  the  buoyancy  in  the  core.  Thus,  this  second  case  appears  a 
much  less  plausible  situation  than  the  first  one. 
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CASE  III:  Balance  between  Convection  and  Conduction  in  the  End 
Region 

In  the  previous  analysis  of  Cases  I  and  II,  situations  were 

considered  in  which  the  driving  force  exists  in  the  core  and  the 

approximate  criteria  were  given  therein  within  which  each  analysis 

2 

is  valid.  Beyond  such  criteria,  i.e.,  for  Ra^A  >>  1,  convection 
becomes  dominant  and  the  diffusion  of  heat  is  negligible  in  the 
core,  thus  most  of  the  temperature  drop  (or  increase)  will  occur  in 
the  end  region  and  will  be  convected  across  the  cavity  to  the 
opposite  end.  This  implies  the  existence  of  thermal  boundary 
layer  in  the  end  region  within  which  most  of  the  temperature  drop 
(or  increase)  occurs.  It  also  implies  the  existence  of  thermal 
boundary  layer  structure  along  the  horizontal  boundaries  in  the 
core  through  which  the  heat  will  be  convected  to  the  opposite  end. 
There  thus  exists  no  driving  force  in  the  core,  instead  it  will  come 
from  the  end  region.  In  what  follows,  consideration  is  given  to 
this  situation  according  to  Pr,  because  physically  for  different 
Pr,  different  force  balances  need  to  be  made. 

3.3.1  Pr  frl  (including  Pr  >>  1) 

3. 3. 1.1  Basic  Equations 

For  the  situation  mentioned  above,  the  heat  transfer  between 
the  two  end  walls  is  dominated  by  convection.  Heat  transfer  by 
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conduction  is  assumed  to  be  important  only  in  thin  boundary  layers 
adjacent  to  the  end  walls.  We  thus  balance  the  convection  and 
conduction  in  the  end  region.  In  addition,  since  all  the  buoyancy 
force  acts  within  the  thermal  boundary  layers  in  the  end  region  and 
the  viscous  effects  are  important  therein  for  Pr  £  1  (including  Pr 
»  1)  because  the  flow  boundary  layers  extend  beyond  the  thermal 
boundary  layers  for  large  Pr,  we  also  balance  the  buoyancy 
and  viscous  forces  in  the  end  region,  by  means  of  which 
the  stretching  parameter,  c  ,  can  be  determined. 


From  Eq.  (2.15),  the  balance  between  convection  and  conduction 
in  the  end  region  can  be  represented  as 
,2 


aL  A 


(3.85) 


From  Eq.  (2.13),  the  balance  between  buoyancy  and  viscous 
forces  in  the  end  region  can  be  represented  as 

6gATi.2H  1  v  L  A? 


*„T  -X 


(3.86) 


Considering  that  the  flow  is  driven  by  the  buoyancy  force  in 
the  end  region,  the  vorticity  would  be  dominant  therein.  It  is  thus 
appropriate  to  represent  the  characteristic  vorticity,  =R,  by 
specifying  the  characteristic  length  2  as  6  ,  as 


4  -v  6. 


(3.87) 
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and  er  ^  (3.88) 

From  (3.85)  n,  (3.87),  we  then  obtain 

vR  %  aRa^  (3.89) 

and  ^  (3.90) 

From  (3.52)  and  (3.90),  we  find 


In  (3.91),  it  is  seen  that  the  end  region  characteristic  length 
6x  is  similar  to  the  familiar  thermal  boundary  layer  thickness  of  a 
vertical  flat  plate. 

Substitution  of  c  and  ;■  into  Eqs.  (2.13)  '*■  (2.15)  then 

K  X  X 

yields 


Considering  the  terms  of  n  and  y  derivatives  in  Eqs.  (3.92) 
(3.94),  the  equations  for  the  end  region  flow  can  be  extracted  as 


(3.95) 


j_  ar*  jj 
* 


2lL 

‘J’7* 


+ 


Jt-  ilL 
*** 


5 


_/  £L 

\ 


tt 


(3.96) 


3(g.  £)  _  Xi.  +.  '  21L 

k^))  ‘  **?  *r 


(3.97) 


With  its  characteristic  length  6x  in  Eq.  (3.91),  the  end  region 
flow  equations  (3.95)  •v  (3.97)  are  similar  to  the  boundary  layer 
equations  of  a  vertical  flat  plate.  Under  similar  boundary  condi¬ 
tions,  the  flow  and  heat  transfer  characteristics  in  the  end  region 
will  thus  be  similar  to  those  in  the  case  of  a  vertical  flat  plate 
except  that  they  are  influenced  by  the  core. 

Considering  the  derivatives  with  respect  to  c  and  y,  the 
equations  which  will  describe  the  core  flow  can  be  written  as 


_L  HI  1 1 

Pr  m.v 


l 


9 (d,  V) 


=  9$  .  1  t ,» ]Q_  'l 

(3.98) 

- - !_  /  +  3L t  \ 

^  +  ar  ) 

(3.99) 

-  .  (  —  /  a*  ^9  .  -3*9  \ 

(3.100) 

Let  us  look  at  the  core  flow  characteristics  based  on  the  di- 
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mensionless  parameters  in  the  Eqs.  (3.98)  ^  (3.100). 


3. 3. 1.2  Core  Flow  Characteristics 


W  <  >• ft2  «  1 

In  Eq.  (3.100),  when  <  1  (or  ARa^  >1)  the  conduction 

H 

terms  are  negligible  compared  to  the  convection  terms.  But  as 
these  terms  are  the  derivatives  of  the  highest  order  in  the  equation, 
for  negligible  horizontal  conduction  term  the  equation  is  sin¬ 
gular  with  respect  to  y  in  the  core.  This  implies  that  there 
exists  a  thin  layer  (it  is  called  here  "horizontal  thermal  layer") 
adjacent  to  the  horizontal  boundaries  within  which  the  vertical 
conduction  term  becomes  important.  This  horizontal  layer  can  be 
estimated  by  the  coordinate  stretching  method.  For  the  consider¬ 
ation  of  the  horizontal  thermal  layer  we  thus  introduce  the  trans¬ 


formation  as 

y  E  ey  y 

where 


(3.101) 

(3.102) 


and  6^  represent  the  horizontal  thermal  layer. 

Substituting  the  derivatives  of  (3.101)  into  Eq.  (3.100)  and 
balancing  the  convection  and  conduction  within  5^,  we  then  have 


1 

ARaHi 


From  (3.102)  and  (3.103),  we  find 


(3.103) 
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or 


<5  'v. 


y 


L 

Ra  i 
H 


(3.104) 

(3.105) 


In  order  that  the  horizontal  thermal  layer  be  distinct. 


6 


/  *  1 

(3.106) 

so  that  AR aH*  >  1 

(3.107) 

Under  the  condition 

(3.107),  to  a  good  approximation, 

the  core 

flow  equation  (3.98)  and 

(3.100)  reduce  to 

J_  _ 

Pr  ' 

90 

(3.108) 

JjjjJtL  - 

0 

(3.109) 

Eqs.  (3.108)  and  (3.109)  will  thus  describe  the  core  configuration 
outside  the  horizontal  thermal  layer,  i.e.,  in  the  mid-core  region, 
under  the  condition  (3.107).  Since  Pr  acts  as  a  parameter  in  Eq. 
(3.108),  we  first  examine  the  core  flow  characteristics  for  high  Pr 
(Pr  >>  1)  and  then  for  moderate  Pr  (Pr  \  l). 

(1)  Pr  »  1 

For  Pr  »  1,  the  inertia  term  in  Eq.  (3.108)  becomes  negligible 
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and  we  have 

||  =  0  (3.110) 

From  (3.110),  two  possible  temperature  profiles  can  be  obtained  as 

6c  =  ec(y)  (3.111) 

or  ec  =  const.  (3.112) 

where  the  subscript  c  is  affixed  to  denote  the  temperature  distri¬ 
bution  outside  the  horizontal  thermal  layer.  Among  the  two  possible 
core  configurations,  it  is  shown  in  Appendix  C  that  the  case  ec  = 
const,  is  inappropriate  to  the  situation  concerned  herein. 

From  (3.109)  and  (3.111),  we  find 

t  =  0  (3.113) 

and  t  ^ (y )  (3.114) 

From  (3.111)  and  (3.114),  it  therefore  can  be  seen  that  for  Pr  >>  1, 
under  the  condition  ARaH*  >  1  there  exists  the  horizontal  thermal 
layer  in  the  core  and  the  temperature  distribution  outside  the 
horizontal  thermal  layer,  i.e.,  in  the  mid-core,  will  be  stratified 
while  the  corresponding  core  flow  structure  will  be  parallel.  This 
parallel  core  flow  structure  was  observed  in  the  experiment  of 
Ostrach  et.  al.  [31]  for  Pr  =  1.38  x  10^,  A  =  0.1  and  RaH  ^  10^> 
but  the  horizontal  walls  are  semi -conducting  therein. 


(1i )  Pr  ^  1 
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When  Pr  •v  1,  (in  fact  Pr  is  about  1  or  slightly  greater 
than  1),  the  inertia  term  in  Eq.  (3.108)  is  not  negligible.  Instead 
the  circulating  flow  will  transport  vorticity  across  the  cavity 
by  the  inertia  and,  as  can  be  conceivable  by  the  singular  behaviour 
in  the  equation,  the  diffusion  of  vorticity  will  be  important  within 
a  layer  (we  call  it  here  "horizontal  viscous  layer"  in  distinction 
to  the  horizontal  thermal  layer)  along  the  horizontal  boundaries. 

This  horizontal  viscous  layer  can  be  estimated  form  the  balance 
between  the  vorticity  transport  and  diffusion  terms  within  that 
layer.  But  in  this  case,  we  have  to  use  a  modified  characteristic 
stream  function  "Vp  instead  of  Vp  ,  because  the  characteristic 
stream  function  has  different  values  according  to  the  different 
characteristic  length  scales  within  which  difference  balances  are 
made.  In  Appendix  D,  the  modified  characteristic  stream  function 
?P  is  estimated.  The  argument  of  the  necessity  of  the  distinction 
between  the  two  characteristic  stream  function  is  also  given  therein. 

For  an  estimate  of  the  horizontal  viscous  layer,  we  introduce 
the  transformation  as 

y  =  ev  y  (3.115) 

where 

6 

ev  (3.116) 

and  6v  is  the  horizontal  viscous  layer. 

Substituting  the  derivatives  of  (3.115)  into  Eq.  (3.98)  with 
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the  modified  characteristic  stream  function  replaced  by  iR,  from 
the  balance  between  the  vorticity  transport  and  diffusion  terms 
within  6  ,  we  have 


1_  ^  vL_ 


*RH  -v 


From  Appendix  0,  as 


vR  a  P  Aa^ 


(3.117) 


(3.118) 


from  (3.116),  (3.117)  and  (3.118),  we  find 

pA 

v  Ran1 


or 


Pr* 


For  distinct  horizontal  viscous  layer. 


so  that  ARau^  >  Pr* 

n 


(3.119) 

(3.12C) 


(3.121) 

(3.122) 


Since  Pr  ^  1,  from  (3.104)  and  (3.119) 

"v  Pr*  %  1  (3.123) 

6y 

Under  condition  (3.122),  the  effect  of  vise  us  shear  in  the 
core  is  confined  to  the  horizontal  viscous  layer,  6y.  Here,  since 
the  flow  is  driven  by  the  boundary  layer  in  the  end  region  and  there 
is  no  other  way  to  induce  any  fluid  motion  in  the  core  except  by 


the  viscous  shear  of  the  end  driven  circulating  flows.  Since 

from  (3.123),  most  of  the  flow  in  the  core  will  thus  circu¬ 
late  through  the  horizontal  viscous  layer  adjacent  to  the  horizontal 
boundaries.  Outside  that  layer,  the  flow  which  may  result  from  the 
entrainment-detrainment  of  the  end  driven  core  circulating  flow 
would  be  of  much  lower  order  than  the  circulating  flow  so  that  the 
motion  therein  would  be  almost  stagnant. 

If  we  compare  the  order  of  vorticity  within  the  horizontal 
viscous  layer  to  that  outside  the  horizontal  viscous  layer,  from 
(3.121)  we  find 


1° 


V"2 

V',2 


■  <r>‘  «  1 


(3.124) 


where  subscripts  i  and  o  denote  the  region  within  and  outside  the 
horizontal  viscous  layer,  respectively. 

Outside  the  horizontal  viscous  layer,  we  may  thus  put 


vc  =  const. 

Then  since  5y  n,  6  ,  from  (3.108)  we 


and  thus  e  =  e  (y) 
c  cKJ ' 


(3.125) 

find 

(3.126) 

(3.127) 


as  indicated  in  Appendix  C. 
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For  =  const.,  there  is  no  convection  at  all  outside 
c  v 

ana  the  energy  equation  (3.1U9)  is  automatically  satisfied. 

From  (3.125)  anc  (3.127),  it  is  thus  seen  that  for  Pr  .  1, 

under  the  condition  ARa^4  >  1  there  exists  the  horizontal  thermal 

layer  in  the  core  and  the  temperature  distribution  in  the  mid-core 

will  be  stratified  while  the  fluid  motion  therein  would  be  almost 

stagnant.  This  stagnant  core  configuration  was  observed  in  the 

experiment  of  Al-Homoud  [32]  for  Pr  =  7.0,  A  =  0.0625  and  Ra^  = 

8  9 

2.0  x  10  ^  2.0  x  10  .  This  shows  good  agreement  with  the  above 
prediction. 

[2]  — -  <  1,  A2  <<-  1 

ARV _ 

Under  this  condition,  the  conduction  term  becomes  important  in 
Eq.  (3.100)  and  from  (3.105),  the  horizontal  thermal  layer 
becomes  of  order 

/  ■v  — — f  /  1  (3.128) 

ARa«- 

There  thus  will  no  longer  exist  distinct  horizontal  thermal  layers 
in  the  core.  Instead  some  horizontal  temperature  gradient  will 
exist  in  the  core  which  may  develop  fluid  motion  in  the  core  in 
addition  to  the  end  driven  circulating  flow.  Further  it  is  supposed 
that  the  thermal  boundary  layer  structure  in  the  end  region  may 
somewhat  be  altered  due  to  the  core  temperature  gradient.  Thus  the 
resulting  flow  driving  mechanism  should  be  modified  from  the  strict 
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end  region  boundary  layer  driven  flow  mechanism.  The  flow  charac¬ 
teristics  in  this  situation  are,  therefore,  supposed  to  lie  between 
those  in  the  flow  regime  of  Case  I  and  in  the  strict  boundary  layer 
driven  flow  regime  of  Case  III.  In  this  sense  this  flow  regime  is 
called  herein  the  "Intermediate  Flow  Regime".  The  condition, 

ARaH*  >  1,  in  (3.107)  would  then  be  a  necessary  condition  for  the 
existence  of  a  distinct  boundary  layer  flow  regime  for  Pr  £  1 
(including  Pr  >>  1).  This  argument  can  be  verified  by  the  experi¬ 
mental  results  of  Imberger  [29],  Al-Homoud  [32]  and  Kamotani  et. 
al  [33]. 

Flow  characteristics  with  the  core  configurations  in  the  inter¬ 
mediate  flow  regime  will  be  considered  in  the  detailed  analysis, 
because  it  is  not  clear  here  whether  its  characteristic  length 
scales  would  be  geometric  length  scales  or  not. 

3.3.2  Pr  <  1  * 

3. 3. 2.1  Basic  Equations 

For  Pr  <  1,  the  flow  boundary  layer  extends  less  than  the 

thermal  boundary  layer  and  the  main  body  of  fluid  can  be  considered 
to  be  inviscid  within  the  thermal  boundary  layer  except  in  the  vi¬ 
cinity  of  end  walls.  Since  all  the  buoyancy  force  acts  within  the 
thermal  boundary  layer,  in  addition  to  the  balance  between  convection 
and  conduction  in  the  end  region,  we  also  balance  the  buoyancy  and 
inertia  forces  in  the  end  region  from  which  the  stretching  parameter, 


*  This  always  includes  the  case  of  Pr  «  1 


ex>  can  be  determined. 

From  Eq.  (2.13),  the  balance  between  buoyancy  and  inertia 
forces  in  the  end  region  can  be  represented  as 


(3.129) 


Then,  with  the  balance  between  convection  and  conduction  in  the 
end  region  in  (3.85),  from  (3.87)  and  (3.129)  we  obtain 


i R  a(PrRaR)i 

and  ex  '■  TP>RaH)i 

From  (3.52)  and  (3.131),  we  find 


(3.130) 

(3.131) 


(3.132) 


Substituting  vD,  e  and  f  into  Eqs.  (2.13)  -v  (2.15),  we  then 

K  X  X 

obtain 


(3.134) 


Considering  the  derivatives  with  respect  to  n  and  y  in  Eqs. 
(3.133)  'v.  (3.135),  the  end  region  flow  equations  can  be  represented 


as 


HULL 


i 


(3.136) 


+  _J _  JiL 


(3.137) 


21111  -  te  .  / 


With  the  characteristic  length  dx  in  Eq.  (3.133),  the  equations 
(3.136)  'v  (3.138)  are  also  similar  to  the  boundary  layer  equations 
of  a  vertical  flat  plate. 

Considering  the  derivatives  with  respect  to  ;  and  y,  the  core 
flow  equations  can  be  represented  as 


2LL12  =  21 

HI,})  'tS 


+ 


Pr 
A  (Pr 


(3.139) 


(3.140) 


21122 

v 


/ 


ACPr^y*  7S 


\ n 


+  JQL  ) 


(3.141) 


I 
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In  the  following,  we  examine  the  core  flow  characteristics 
based  on  the  dimensionless  parameters  among  the  core  flow  equations. 


3. 3. 2. 2  Core  Flow  Characteristics 


fl]  - -  <  i,  A2  «  1 

A(PrRaH)4 


A(PrRaH V 


<  1,  both  the  viscous  and  heat  diffusion  terms 


become  negligible  in  Eqs.  (3.139)  and  (3.144)  and  as  in  the  previ¬ 
ous  analysis,  the  horizontal  layers  can  be  estimated  by  the  coordin¬ 
ate  stretching  method. 

Substituting  the  derivatives  of  (3.101)  into  Eq.  (3.141)  and 
balancing  the  convection  and  conduction  within  6  ,  we  have 


(3.142) 


From  (3.102)  and  (3.142),  we  then  find 


(3.143) 


In  order  that  the  horizontal  thermal  layer  be  distinct, 


(3.144) 


or  A(PrRaH)i  >  1 


(3.145) 
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Under  the  condition  (3.145),  to  a  good  approximation,  outside  the 
horizontal  thermal  layer  6^,  i.e.,  in  the  mid-core,  the  core  flow 
equaitons  (3.139)  and  (3.141)  reduce  to 


2LL12 

V?,») 


(3.146) 


ilMi  - 
'&(*,**) 


o 


(3.147) 


For  Pr  <  1,  as  can  be  seen  from  the  works  of  Ostrach  [43]  and 
Sparrow  and  Gregg  [44],  the  thickness  of  velocity  boundary  layers 
along  the  end  walls  is  less  than  (or  about  equal  to)  that  of  ther¬ 
mal  boundary  layers.  In  the  core,  it  is  also  expected  that  the 
horizontal  viscous  layer  is  less  than  (or  about  equal  to)  the 
horizontal  thermal  layer.  In  this  situation,  as  mentioned  in 
Appendix  D,  the  value  of  characteristic  stream  function  within  the 
horizontal  thermal  layer  will  be  identical  to  that  within  the 
horizontal  viscous  layer.  There  thus  is  no  need  to  modify  the 
characteristic  stream  function  for  the  estimation  of  horizontal 
viscous  layer.  Then,  by  introducing  the  transformation  of  (3.115) 
into  Eq.  (3.139)  and  from  the  balance  between  the  vorticity  trans¬ 
port  and  diffusion  terms  within  the  horizontal  viscous  layer  £  ,  we 
find 
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For  distinct  horizontal  viscous  layer 

^  <  1  (3.149) 

so  that  A{PrRau)i  >  Pr  (3.150) 

n 

From  (3.144)  and  (3.148),  it  can  be  shown  that 

6  <  6  (3.151) 

v  y 

for  Pr  <  1. 

Under  condition  (3.150),  the  effect  of  viscous  shear  in  the 
core  will  be  confined  to  the  horizontal  viscous  layer  and  since 
6^  <  6^  in  (3.151),  most  of  the  flow  in  the  core  will  thus  circulate 
through  the  horizontal  thermal  layer.  Outside  the  horizontal  ther¬ 
mal  layer,  as  mentioned  earlier  in  the  previous  analysis,  the 
fluid  motion  would  remain  nearly  stagnant.  We  may  thus  put 

ic  =  const.  (3.152) 

in  the  mid-core  region. 

Then,  from  Eq.  (3.146)  we  find 

||  =  0  (3.153) 

and  thus  ec  =  ec  (y)  (3.154) 

as  indicated  in  Appendix  C. 

It  is  thus  seen  from  (3.152)  and  (3.154)  that  for  Pr  <  1, 
under  the  condition  A(PrRa„)*  >  1,  with  the  distinct  horizontal 
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thermal  layer  in  the  core,  the  temperature  distribution  will  be 
stratified  in  the  mid-core  while  the  fluid  motion  therein  would 
remain  almost  stagnant.  Unfortunately  no  experiments  have,  as  yet, 
been  available  for  Pr  <  1  to  verify  the  above  prediction. 


[2] 


- - r  l  1,  A2  «  1 

A(PrRaH)* 


Under  this  condition,  from  Eq.  (3.144),  conduction  becomes 
important  in  the  core  and  from  Eq.  (3.146),  the  horizontal  thermal 
layer,  6  ,  becomes  of  order 


H  A(PrRaH)i 


> 


(3.155) 


Thus,  there  will  no  longer  exist  the  distinct  horizontal  thermal 
layer  in  the  core,  instead  some  core  fluid  motion  can  be  induced  by 
the  horizontal  temperature  gradients  in  the  core  in  addition  to  the 
boundary  layer  driven  circulating  flows.  Also  the  thermal  boundary 
layer  structure  in  the  end  region  and  thus  the  resulting  flow  dri¬ 
ving  mechanism  will  be  modified  according  to  the  core  temperature 
gradients.  This  flow  regime  also  corresponds  to  the  "intermediate 
flow  regime"  mentioned  previously  and  will  be  considered  later  in 
the  detailed  analysis.  The  condition  (3.147),  i.e.,  A(PrRaH)*  >  1, 
would  thus  be  the  necessary  condition  for  the  distinct  boundary 
layer  flow  regime. 

In  this  chapter,  we  have  studied  the  global  flow  characteris¬ 
tics  according  to  three  different  balances:  the  balance  between 


buoyancy  and  viscous  forces  in  the  core  in  Case  I,  the  balance  be¬ 
tween  buoyancy  and  inertia  forces  in  the  core  in  Case  II  and  the 
balance  between  convection  and  conduction  in  the  end  with  the  addi¬ 
tional  force  balance  according  to  Pr.  In  what  follows,  considera¬ 
tions  will  be  given  both  to  the  core  flow  structure  of  the  inter¬ 
mediate  flow  regime  and  to  the  possible  flow  subregimes  in  the 
boundary  layer  flow  regime  of  Case  III. 


CHAPTER  IV 


DETAILED  VIEW  OF  THE  FLOW  CHARACTERISTICS 

In  the  previous  analysis,  the  core  flow  characteristics  were 
studied  globally  according  to  three  different  cases.  In  Case  I, 
the  balance  was  made  between  buoyancy  and  viscous  forces  in  the  core 
and  the  core  flow  pattern  was  shown  to  be  parallel  whereas  the  core 
flow  pattern  was  shown  to  be  non-parallel  in  Case  II,  in  which  the 
balance  was  made  between  buoyancy  and  inertia  forces  in  the  core. 

In  Case  III,  the  balance  was  made  between  convection  and  con¬ 
duction  in  the  end  region.  In  addition,  a  balance  was  made  in  the 
end  region  between  buoyancy  and  viscous  forces  for  Pr  1  (including 
Pr  >>  1)  and  between  buoyancy  and  inertia  forces  for  Pr 1.  In  Case 
III,  when  there  exist  distinct  horizontal  thermal  layers  in  the 
core  the  flow  was  shown  to  be  driven  by  the  boundary  layers  in  the 
end  region  (called  boundary- layer  flow  regime).  In  this  situation 
the  core  flow  structure  was  shown  to  be  parallel  for  Pr  >>  1,  while 
the  fluid  motion  is  expected  to  remain  almost  stagnant  in  the  mid¬ 
core  both  for  Pr  ^  1  and  Pr  <  1.  In  addition,  when  the  flow  is 
driven  by  the  end,  there  may  appear  some  flow  subregimes  in  the  core 
due  to  the  local  interactions  of  the  end  and  horizontal  boundary 
layers.  The  occurrence  of  such  flow  subregimes  will  change  the 
global  core  flow  patterns,  and  their  characteristic  length  scales 
would  also  be  different  from  the  geometric  length  scales.  It  thus 
becomes  important  to  be  able  to  predict  the  conditions  under  which 
flow  subregimes  are  present  and  to  estimate  their  length  scales  to 
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determine  the  core  flow  pattern. 

In  the  analysis  of  Case  III,  in  addition  to  the  boundary  layer 
flow  regime  the  term  "Intermediate  Flow  Regime"  was  defined  in  the 
sense  that  its  flow  and  heat  transfer  characteristics  would  lie 
between  the  flow  regime  wherein  the  flow  is  driven  in  the  core 
[Cases  I  and  II]  and  that  the  driver  is  in  the  end  region  [boundary 
layer  flow  regime  in  Case  III].  For  such  a  flow  regime,  since  it 
was  not  clear  whether  the  geometric  length  scales  are  appropriate 
for  the  description  of  core  flow  structure,  no  global  analysis  of 
the  flow  pattern  was  made. 

In  what  follows,  consideration  will  be  given  to  the  prediction 
of  the  possible  flow  subregimes  in  the  boundary  layer  flow  regime 
and  also  to  the  prediction  of  the  core  flow  pattern  in  the  inter¬ 
mediate  flow  regime,  both  of  which  could  not  be  predicted  from  the 
global  analysis. 

4.1  Prediction  of  Flow  Subregimes  in  the  Boundary  Layer  Flow 
Regime 

4.1.1  Working  Equations 

The  characteristic  length  scales  of  either  the  core  flow 
structure  in  the  intermediate  flow  regime  or  the  expected  flow  sub¬ 
regimes  in  the  boundary  layer  flow  regime  may  be  different  from  the 
geometric  length  scales.  We  thus  replace  the  geometric  length 
scales  hc  and  in  the  dimensionless  equations  (2.13)  and  (2.15), 
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which  are  to  be  specified  in  the  course  of  analysis.  In  addition, 
as  it  is  seen  from  the  global  analysis  that  the  basic  equations 
show  singular  behaviour  at  the  horizontal  boundaries,  in  order  to 
account  for  that  we  introduce  other  multiple  scales  for  the  vertical 
coordinate  y  as 


y  =  y. 


=  JL 


(4.1) 


where 


ch  %  h 


(4.2) 


and  6h  is  an  arbitrary  horizontal  layer  to  be  properly  specified  in 
the  analysis.  The  characteristic  stream  function  is  also  to  be 
specified  according  to  6h  concerned. 

Introducing  the  derivatives  of  multiple  scales  (4.1)  into  Eqs. 
(2.13)  and  (2.15)  and  replacing  H  and  L  by  hc  and  *  ,  we  then  have 


4-4.  'J-  t  ^ 


1LL1L  ,  >  m.p  ,  Vt.t)  - 

MV.V  4  9(5.?)  mt)  w.  9^ 


c 


i  J2JL- 


(4.3) 


_L_L  HM2. 


+  J-2LAt± 

u  m.}) 


4,  + 


-L  ^Xr  +  jLfi-'i 


/V  re. 

U*  *7* 

(4.4) 
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h 

where  A  =  and  6  replaces  t  in  Eq.  (4.3) 
c  'c  x 

4.1.2  Flow  Subregimes  near  the  End  Region 

It  was  mentioned  earlier  that  flow  subregimes  are  expected  in 
the  boundary-layer  flow  regime  due  to  the  interaction  of  the  end  and 
horizontal  boundary  layers.  We  first  consider  the  situation  in 
which  the  end  region  boundary-layer  flow  turns  around  the  upper  cor¬ 
ner  of  the  hot  wall  (or  lower  corner  of  the  cold  wall).  In  this 
situation,  since  the  velocity  has  different  values  near  the  two 
rigid  boundaries  at  the  corner,  the  gradients  of  velocity  become 
very  large  so  that  the  viscous  friction  is  supposed  to  be  very 
large  at  the  horizontal  boundaries  near  the  corner.  In  this  sense, 
the  present  situation  is  physically  similar  to  the  case  of  two- 
dimensional  flow  in  a  corner  due  to  one  rigid  plane  sliding  on 
another  with  constant  inclination  angle  ,  although  the  flow 
herein  is  driven  by  the  buoyancy  force  within  the  thermal  boundary 
layer  at  the  end  walls. 

When  the  heated  boundary-layer  flow  turns  around  the  corner, 
the  streamlines  of  the  heated  flow  can  be  deflected  due  to  the 
large  viscous  friction  there.  If  the  streamline  deflection  is 
large,  some  of  the  unheated  flow  entrained  by  the  viscous  shear  of 
heated  boundary-layer  flow  can  be  blocked  so  that  it  does  not  reach 
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the  opposite  end.  In  this  way,  closed  streamlines  could  apoear  as 
flow  subregimes  such  as  secondary  cells.  However,  this  does  not 
mean  that  the  heated  boundary-layer  flow  is  also  closed.  Even 
though  the  heated  flow  is  deflected  near  the  corner,  it  will  soon 
rise  up  and  will  circulate  along  the  horizontal  boundaries  develop¬ 
ing  the  horizontal  thermal  layers  in  the  core.  Otherwise,  there 
will  not  exist  distinct  horizontal  thermal  layers  in  the  core,  so 
that  the  thermal  boundary  layer  structure  at  the  end  walls  will  be 
altered  and  the  distinct  boundary-layer  flow  regmie  will  not  exist. 
As  a  matter  of  fact,  secondary  cells  near  the  end  region  were  ob¬ 
served  in  some  experiments  and  our  attention  here  is  focused  on  the 
prediction  of  such  secondary  cells  considering  the  proper  force 
balances.  The  inferred  flow  pattern  for  the  above  situation,  which 
is  actually  based  on  the  experimental  observations,  is  approximately 
as  indicated  in  Fig.  2. 

If  we  consider  the  globally  clear  driving  mechanism  of  the 
present  situation,  it  can  be  inferred  that  the  extent  of  the  stream¬ 
line  deflection  will  largely  depend  on  the  interaction  of  the  iner¬ 
tia  of  the  buoyancy  driven  boundary-layer  flow  and  the  viscous 
friction  within  an  arbitrary  horizontal  layer  6,  along  the  horizon¬ 
tal  boundaries.  Since  the  heated  boundary  layer  flow  will  circulate 
across  the  top  of  the  enclosure  through  the  horizontal  thermal  layer 
6^,  specifying  6 as  6^  in  (4.2),  from  Eq.  (4.3),  the  proper  force 
balance  of  the  boundary  layer  inertia  and  horizontal  viscous  fric¬ 
tion  for  that  interaction  can  be  represented  as 


re  2. 


Approximate  Representati  >n  of  the  Flow 
Pattern  of  Concern  near  the  End  Region 
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VJt, 


1 


S* 


and  this  gives 


,  e»  Vc 

£  «v  -  - 

c  s  v 


(4.5) 


(4.6) 


Since  the  characteristic  values  of  vD  and  c  in  (4.6)  differ 

K  X 

according  to  Prandtl  numbers  in  the  boundary-layer  flow  regime, we 
first  consider  the  case  of  Pr  ^  1  (including  Pr  >>  1)  and  then  that 
of  Pr  <  1  (as  mentioned  earlier,  this  always  includes  the  case  of 
Pr  «  1). 


(i )  Pr  1  (including  Pr  >>  1 ) 


In  this  case,  substituting  the  characteristic  values  vR 

aRa*  and  ev  'v  given  in  (3.89)  and  (3.90)  into  (4.6),  the 

M  x  RaH* 

result  is 


Ra, 


*  n,  e, 

c  h 


PrA 


(4.7) 


To  determine  the  length  scales  hc  and  a  in  (4.7)  explicitly, 
one  more  relation  is  required.  Considering  the  globally  clear 
physical  situation,  however,  no  additional  balance  appropriate  to 
the  present  situation  seems  to  exist.  But,  based  on  the  experi¬ 
mentally  observed  secondary  cell  patterns  as  shown  in  Fig.  2,  it 
seems  reasonable  to  suppose 
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hc  %  H 

From  (4.7),  we  then  obtain 


Ra, 


Pr 


L  = 


A 


Gr, 


Lh  /  Pr 


(4.8) 


(4.9) 


For  large  streamline  deflection  of  the  heated  boundary  layer 
flow  near  the  corner,  the  horizontal  thermal  layer  6y  can  be  re¬ 
garded  as  of  order 

6y  *  0(H) 

6  & 

or  Eh  "  /  *  "  1 

c 

From  (4.9)  and  (4.11),  we  then  have 


(4.10) 

(4.11) 


In  (4.12),  when  the  arbitrary  horizontal  length  scale  is 
less  than  L,  i.e.,  when  —  <  1,  flow  subregimes  are  expected  to 
appear.  This  statement,  though,  is  not  the  necessary  and  sufficient 
conditions  for  the  occurrence  of  the  flow  subregimes,  but  it  at 
least  gives  the  necessary  condition. 

From  (4.12),  for  <  1 

Pr*  >  GrH* 


(4.13a) 
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or  Pr  >  RaR*  (4.13b) 

For  large  values  of  RaR,  the  condition  (4.13)  can  thus  be  satisfied 
for  Pr  >>  1. 

The  types  of  flow  subregimes  will  depend  on  the  length  scale 
t  .  From  (4.8)  and  (4.12),  it  is  seen  that 


He  %  /Pr_  A 
*c  J  GrH 

(4.14) 

When 

APr‘  %  Gr^ 

(4.15a) 

or 

APr  Ra  * 

H 

(4.15b) 

we  find 

r  *  1 

(4.16) 

This  indicates  that  the  flow  subregime  will  appear  as  secondary  cells 
near  the  end  region.  Considering  the  global  flow  characteristics  in 
the  core  for  Pr  >>  1,  however,  patterns  of  multi-cells  are  not  ex¬ 
pected  to  the  present  situation. 

( i  i )  Pr  <  1 

If  we  substitute  the  characteristic  values  of  vR  %  A(PrRaH)* 
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and  ex  ^  '['p'rR'a  ~JT  91ven  in  (3.130)  and  (3.131)  into  (4.6)  and 
H 

also  supposing  hc  %  H,  we  obtain 

‘c  -  ch2  GrH*  L  (4.17) 

For  large  streamline  deflection,  from  (4.10)  and  (4.11)  we  may 
also  put  'v  1.  From  (4.17),  we  then  have 

£c  %  GrH*  L  (4.18) 

Since  Gr^  >>  1  for  Pr  <  1  in  the  boundary  layer  flow  regime,  from 
(4.18)  we  find 


(4.19) 


This  indicates  that  no  flow  subregimes  are  expected  near  the  end 
region  for  Pr  <  1. 

From  the  above  analysis,  flow  subregimes  near  the  end  region 
are  expected  only  for  Pr  >>  1.  This  is  due  to  the  fact  that  in  the 
case  of  Pr  •v  1  and  Pr  <  1,  the  inertia  of  the  boundary  layer  flow 
is  strong  enough  to  overcome  the  viscous  friction  near  the  corner 
while  it  turns  around  there,  so  that  no  appreciable  streamline 
deflection  occurs  near  the  end  region. 

Experimentally  secondary  cells  near  the  end  region  were  ob- 

2 

served  in  the  work  of  Ostrach  et.  al .  [31]  for  A  =  0.2,  (RaHA)  n, 
10*3  and  Pr  =  1.38  x  103  and  Kamotani  et.  al  [33]  for  A  =  0.2, 

Ra^  =  4.8  x  10^  and  Pr  =  965.  In  [31],  secondary  cells  appeared 
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for  somewhat  higher  values  of  Ra^  than  the  predicted  value.  But  in 
that  case,  the  horizontal  boundaries  are  semi -conducting  rather  than 
adiabatic.  The  results  in  [33]  show  very  good  agreement  with  the 
above  prediction.  In  the  experiment  of  Simpkins  and  Dudderar  [38], 
secondary  cells  were  also  observed  and  the  critical  Rayleigh  number 
for  the  occurrence  of  secondary  cells  is  given  by 

(RaH)c  A3/4  %  6.4  x  105  ±  10%  (4.20) 

which  has  an  aspect  ratio  dependence  for  the  range  of  0.1  *  A  £  4. 
The  condition  (4.20)  does  not  show  good  agreement  with  the  results 
of  [31]  and  [33].  In  the  experiments  of  [31],  [33]  and  [38],  when 
secondary  cells  appeared,  there  was  only  one  in  each  end  of  the 
cavity  and  the  core  flow  was  not  parallel.  No  experiments  have 
been  conducted,  as  yet,  for  A  £  0.1  with  adiabatic  horizontal  boun¬ 
daries  under  the  condition  indicated  in  (4.15)  for  the  existance  of 
secondary  cells.  Since  for  Pr  >>  1,  the  global  analysis  for  the 

parallel  core  flow  structure  in  the  boundary  layer  flow  regime 

2 

yields  only  the  condition  A  <<  1,  it  would  be  interesting  to  make 
observation  under  the  condition  (4.15). 

In  the  numerical  work  of  Jhaveri  and  Rosenberger  [40],  for 
A  =  0.1,  Sc  =  0.5,  Pr  =  0.7  and  Gr^  =  3.3  x  104,  the  resulting  flow 
profile  shows  the  onset  of  a  recirculation  cell  at  each  interface. 

It  should  be  noted  that  this  cell  is  due  to  the  coupling  of  mass 
flux  at  the  interface  with  the  transport  quantities. 
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4.1.3  Flow  Subregimes  in  the  Core 

In  the  preceding  analysis,  secondary  cells  were  shown  to  be 
possible  near  the  end  region  for  Pr  »  1.  Here  the  possibility  of 
flow  subregimes  is  going  to  be  examined  away  from  the  end  region, 
i.e.,  in  the  core  by  considering  the  interaction  of  the  inertia 
and  viscous  forces  of  the  end  driven  core  circulating  flow  along 
the  long  horizontal  boundaries.  This  interaction  was  considered  in 
the  previous  global  analysis  for  Pr  ^  1.  From  that  consideration, 
it  was  shown  that  the  horizontal  viscous  layer  is  distinct  in  the 
core  under  the  condition  ARa^*  >  Pr*  and  the  corresponding  flow 
pattern  remains  nearly  stagnant  in  the  mid-core.  Since  the  con¬ 
dition  required  is  ARa^*  >  1  for  the  existence  of  a  distinct  hori¬ 
zontal  thermal  layer  in  the  core,  the  question  naturally  arises 
what  the  core  flow  pattern  would  be  like  under  the  condition 

1  <  ARa*  <  Pr*,  wherein  the  horizontal  thermal  layer  is  distinct 
H 

but  the  horizontal  viscous  layer  is  of  order,  6y>H.  The  present 
analysis  is  thus  concerned  in  predicting  the  core  flow  pattern  in 
that  situation. 

In  the  analysis  the  possibility  is  implicitly  disregarded  for 
Pr  >>  1,  because  in  that  case,  the  inertia  is  negligible  and  glob¬ 
ally  the  core  flow  structure  was  shown  to  be  parallel  unless  there 
appear  secondary  cells  near  the  end  region. 

If  the  balance  is  made  for  the  interaction  between  the  inertia 
and  viscous  vorces  within  the  arbitrary  horizontal  layer  6h,  from 
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Eq.  (4.3)  it  can  be  represented  as 


eh 


»Rhc 


(4.21) 


Since  we  are  considering  the  momentum  balance  of  the  circula¬ 
ting  flow,  the  arbitrary  horizontal  layer  6^  can  be  regarded  as  the 
horizontal  viscous  layer  6y.  The  characteristic  stream  function 
Tn  in  (4.21)  thus  needs  to  be  replaced  by  the  modified  stream  func¬ 
tion  v-  for  6  .  We  first  consider  the  cis>  of  Pr  ■.  1  and  then  that 

K  V 

of  Pr  <  1. 


(i)  Pr  ^  1 


Replacing  by  yr  in  (4.21)  and  substituting  the  characteris¬ 


tic  value  of  yR  from  Appendix  D,  we  find 


Ra, 


*c  %  eh 


PrJ 


(4.22) 


Here  we  also  need  one  more  relation  to  determine  the  length 

scales  h  and  j i  explicitly.  However,  there  does  not  seem  to  exist 
c  c 

an  appropriate  balance  for  the  present  situation.  As  was  done  pre¬ 
viously,  if  we  suppose  hc  *v  H,  for  large  streamline  deflection,  as 
in  (4.11)  we  may  put 

6.. 

(4.23) 


From  (4.22)  and  (4.23),  we  then  obtain 
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Ra 


H 


Pr 


i 


(4.24) 


For  -j— -  <  1 ,  from  (4.24) 
i 


Ra 


H 


Pr 


i 


A  <  1 


(4.25c) 


or 


ARa, 


<  PrJ 


(4.25b) 


Since  ARaH  >  1  for  the  boundary  flow  regime,  from  (4.25)  flow 
subregimes  are  expected  in  the  core  under  the  condition 


RV 


<  A 


Pr5 


Ra, 


(4.26) 


when  Pr  %  1,  from  (3.123)  «5y  *  <y>  and  thus  for 

{v  %  H  in  (4.23),  is  also  of  order  0(H).  In  this  situation,  no 
distinct  horizontal  thermal  layer  exists  in  the  core  and  the  flow 
regime  corresponds  to  the  intermediate  flow  regime.  Therefore  in 
the  boundary  layer  flow  regime,  from  (4.26),  as  long  as  ARa*  >  1, 
no  flow  subregimes  are  expected  in  the  core  for  Pr  £  l,  but  they  are 
expected  for  Pr  slightly  greater  than  1  (for  example,  water). 

The  types  of  flow  subregimes,  as  mentioned  earlier,  will  depend 
on  the  length  scale  i^.  From  (4.25),  we  find 


Pr  * 


ARa, 


^  n 


(4.27) 
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and 


h 

I 


£ 

C 


Pr3 


Ra, 


(4.28) 


where  n  denotes  an  arbitrary  integer  approximating  the  ratio  of 

lc 

From  (4.27),  the  value  of  n  will  indicate  the  patterns  of  flow 

subregimes  in  the  core.  For  example,  when  n  =  3  this  may  imply  the 

appearance  of  three  cells  in  the  core.  However,  since  ARa^  >  1 

and  the  case  of  Pr  >>  1  is  neglected  here,  the  ratio  of  —  in  (4.27) 

c 

cannot  be  a  much  larger  value  than  1.  Further,  from  (4.28),  for 
large  values  of  Ra^  we  find 


H 


<  1 


(4.29) 


for  Pr  v  1.  Therefore,  patterns  of  multi -cel Is  are  not  expected 
to  appear  in  this  situation. 


(ii)  Pr  <  1 

In  this  case,  from  (3.151)  it  is  seen  that  the  horizontal  ther¬ 
mal  layer  extends  beyond  the  horizontal  viscous  layer,  i.e.,  > 

6y .  Thus  for  6y  -v  H  in  (4.23),  <5^  is  of  order,  >  H.  Then 
there  does  not  exist  a  distinct  horizontal  thermal  layer  in  the 
core  and  the  flow  regime  in  this  situation  corresponds  to  the  in¬ 
termediate  flow  regime.  For  Pr  <  1,  no  flow  subregime  are,  there¬ 
fore,  expected  in  the  core  as  well  as  near  the  ^nd  region  in  the 
boundary  layer  flow  regime. 
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Experimentally  a  single  slow  central  circulating  flow  subregime 
in  the  core  was  observed  in  the  work  of  Imberger  [29]  for  Pr  =  7, 

O 

Rau  =  1.11  x  10  and  A  =  0.02.  This  shows  good  agreement  with  the 
n 

above  prediction.  No  experiments  with  the  fluids  of  Pr  <  1  have 
been  conducted  in  the  boundary  layer  flow  regime,  as  yet,  for  the 
present  configuration. 

Summarizing  the  above  analysis,  in  the  boundary  layer  flow 
regime,  flow  subregimes  such  as  secondary  cells  are  possible  near 
the  end  region  for  Pr  >>  1  whereas  they  are  possible  in  the  core  for 
Pr  -u  1  (except  Pr  >  1).  No  *iow  subregimes,  however,  are  expected 

for  fl ui ds  wi th  Pr  <  1 . 

4.2  Prediction  of  Core  Flow  Pattern  -in  the  Intermediate  Flow  Regime 
4.2.1  Working  Equations 

In  the  intermediate  flow  regime,  there  does  not  exist  distinct 
horizontal  thermal  layer  in  the  core.  Thus-  there  would  be  no  sin¬ 
gular  behaviour  of  the  basic  equations  with  respect  to  thn  horizon¬ 
tal  boundaries  so  that  it  is  not  needed  to  introduce  the  other 
multiple  scales  with  respect  to  y  as  in  (4.1).  In  this  situation 
the  working  equations  (4.3)  and  (4.4)  simply  reduces  to 


L'iil.f).  +  20L12 


(4.30) 
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J_  'MAI  +  we.*) 


i-  3. 


A.  .XL- 

Ck  973? 


+  a*  XL.  a- XL) 
‘ jr  sy J 


(4.31) 


4.2.2  Prediction  of  Core  Flow  Pattern 

The  balances  in  Eqs.  (4.30)  and  (4.31)  which  will  properly 
describe  the  fluid  physics  in  the  intermediate  flow  regime  are  not 
readily  apparent.  However,  as  mentioned  previously,  it  is  supposed 
that  the  flow  and  heat  transfer  characteristics  in  the  intermediate 
flow  regime  would  lie  between  the  flow  regime  wherein  the  flow  is 
driven  in  the  core  and  that  wherein  the  flow  is  driven  by  the  boun¬ 
dary  layer  in  the  end  region.  In  othei  nrds,  it  may  be  said  that 
the  balances  made  in  both  flow  regimes  can  now  be  made  in  the  inter¬ 
mediate  flow  regime.  It  is  thus  from  this  arguement  that  the  balan¬ 
ces  are  going  to  be  made  in  the  following  analysis. 

In  the  flow  regime  wherein  the  flow  is  driven  in  the  core, 
the  force  balances  are  made  in  the  core  either  between  the  buoyancy 
and  viscous  forces  or  between  the  buoyancy  and  inertia  forces 
according  to  Pr.  In  the  end  driven  flow  regime,  the  balances  are 
made  in  the  end  region  either  between  the  buoyancy  and 

viscous  forces  or  between  the  buoyancy  and  inertia  forces.  In 
addition,  in  the  former  case  the  heat  convection  and  conduction 


1 
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terms  become  of  equal  order  in  the  core  for  higher  value  of  Rah 
(i.e.,  for  Ra^A  -v  1  in  Case  1)  whereas  both  terms  are  balanced 
in  the  end  region  in  the  latter  case.  Therefore,  basically  two 
balances  can  be  made  herein  in  both  the  end  and  core  regions  either 
from  the  force  balances  or  from  the  heat  balances.  In  the  following, 
we  first  make  the  heat  balances  in  both  end  and  core  regions. 

Based  on  the  heat  balances,  the  force  balances  will  then  be  made 
appropriately  and  we  will  see  what  the  results  imply.  Other  force 
balances  will  also  be  considered  if  necessary. 

From  Eq.  (4.31),  the  balance  between  convection  and  conduction 
in  the  end  region  can  be  represented  as 


and  in  the  core 


1  -v 


'Rhc 


(4.33) 


Since 


from  (4.32)  and  (4.33),  we  obtain 


(4.34) 


(4.35) 


and 


(4.36) 
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The  explicit  forms  of  the  characteristic  quantities  in  (4.35) 
and  (4.36)  can  be  determined  from  the  length  scales  of  hc  and  »'c 
obtained  by  proper  force  balances  for  the  motions  induced  by  both 
the  core  and  end  temperature  gradients.  From  the  global  analysis, 
it  is  evident  that  the  Pr  acts  as  a  parameter.  Thus  this  will  be 
done  according  to  the  Pr. 

(i)  Pr  1 


In  this  case,  the  balance  between  the  buoyancy  and  viscous 
forces  are  appropriate.  We  first  consider  the  balance  in  the  end 
region.  From  Eq.  (4.30),  this  balance  can  be  represented  as 


.llfTi/A...  A  _  A* 


(4.37) 


From  (4,35)  ,(4.36)  and  (4.37),  we  find 


( 4 . 38 a  ) 


or 


(4 - 


(4.38b) 


In  order  to  determine  the  length  scales  hc  and  ;  explicitly, 
one  more  force  balance  is  required.  Firstly,  if  we  take  the  balance 
between  the  buoyancy  and  viscous  forces  in  the  core,  from  Eq.  (4.30) 
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this  can  be  represented  as 

jglfl&fc  ^  jlL. 


(4.39) 


From  (4.35),  (4.36)  and  (4.39),  we  then  have 


~  if w  <4-40) 

which  is  identical  to  the  relation  (4.38).  Thus  this  force  balance 
does  not  give  any  more  information. 

Secondly,  if  we  consider  the  possible  interaction  of  the  vis¬ 
cous  forces  between  the  end  and  core  region,  one  additional  force 
balance  can  be  inferred,  that  is,  the  balance  between  the  viscous 
force  in  the  end  and  the  viscous  force  in  the  core.  From  Eq.  (4.30) 
this  balance  can  be  represented  as 


This  gives 

£x  '  1 

or  «„  -  «£ 

From  (4.36)  and  (4.43),  we  find 


(4.41) 


(4.42) 

(4.43) 


(4.44) 


so  that  *  ^  h„ 

c  c 


(4.45) 


From  (4.38)  and  (4.45),  we  then  obtain 


u  3  av _ 

hc  egAt 


and  this  shows 
h. 


(4.46) 


H 


Ra 


1/3 


(4.47) 


H 


Since  A^Ra^  »  1  for  the  intermediate  flow  regime,  for  «  1, 
Ra^  is  a  large  value.  From  (4.47),  we  thus  find 


h£<  1 


(4.48) 


The  relat.ons  (4.45)  and  (4.48)  indicate  the  existence  of  very 

small  cells  in  the  core.  No  observation  of  such  cells  has,  as  yet, 

been  reported  within  the  parametric  ranr^s  of  concern  here  (i.e., 

2  ^ 

A  Ra„  »  1,  ARaH  £  1),  although  the  data  are  not  extensive. 
Physically,  the  appearance  of  very  small  cells  in  the  core,  in 
general,  does  not  seem  to  be  a  plausible  situation,  and  thus  the 
balance  between  the  viscous  forces  in  the  end  and  core  region  seems 
to  be  inappropriate  for  the  present  situation. 

If,  as  previously,  we  suppose 


h„  H 
c 


(4.49) 


from  (4.38),  we  obtain 
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*c  *  RaH  H 


(4.50) 


From  (4.49)  and  (4.50),  Vp  and  6x  in  (4.35)  and  (4.36)  can  then  be 


represented  as 


*R  '  “RV 


%  -“-r 

RaH 


(4.51) 


(4.52) 


From  (4.51)  and  (4.52),  it  is  seen  that  the  expression  for  vR 
and  6x  herein  are  similar  to  those  of  the  boundary-layer  flow  re¬ 
gime.  This  may  imply  that  the  flow  characteristics  in  the  inter¬ 
mediate  flow  regime  concerned  here  are  rather  closer  to  those  in 
the  boundary-layer  flow  regime.  However,  since  A RaRi  £  1  in  the 
intermediate  flow  regime,  from  (4.50)  we  find 


^  ARa^  ^  1 


(4.53) 


and  for  ARa^  <  1 


<  1, 


(4.54) 


which  indicates  the  existence  of  flow  subregimes  in  the  core  due  to 
the  core  temperature  gradients.  The  length  scale  «c,  which  is  a 
certain  fraction  of  L  in  (4.54),  may  thus  indicate  what  type  of  flow 
subregimes  will  appear. 

When  ARaH*  <  1,  from  (4.50)  we  may  write 


< 
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L 

t 


% 


1 

ARaH* 


% 


n 


(4.55) 


where  n  is  an  arbitrary  integer  approximating  the  ratio  of  y-  . 

2  2  c 

For  large  values  of  Ra^  such  that  Ra^A  »  1  for  A  «  1, 
from  (4.50)  we  also  find 


h 

I 


c 

c 


< 


1 


From  (4.55)  and  (4.56),  since 


(4.56) 


'L 


£  An 


(4.57) 


we  then  find 

1  <  n  <  j  (4.58) 

This  indicates  that  patterns  of  multi-cells  are  possible  in  the  core 

2 

for  small  values  of  A,  i.e.,  for  A  <<  1. 

*  t- 

When  ARa^  ^  1,  from  (4.50)  ^  1.  In  this  situation,  no 

apparent  flow  subregimes  are  expected  and  it  is  not  readily  clear 

what  type  of  flow  patterns  will  appear  in  the  core, 
i  lc 

When  ARa^  >  1,  —  >  1  and  thus,  no  flow  subregimes  are  expected 
This  is  the  case  of  boundary- layer  flow  regime.  It  is  evident  from 
the  global  analysis  that  no  flow  subregimes  are  expected  by  the  core 
temperature  gradient  in  the  boundary-layer  flow  regime. 
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Few  experimental  data  available  for  the  present  situation  are 
found  In  the  Imberger's  experiment  [29].  In  [29],  when  ARaH*  ^  1, 
e.g.,  Pr  =  7.0,  RaH  =  8.04  x  10®  and  A  =  0.02,  some  part  of  the  core 
flow  turns  around  before  it  reaches  the  end  region.  For  ARaH*  <  1, 
however,  the  flow  structure  is  parallel  in  the  core  and  no  flow 
subregimes  are  observed.  Nevertheless,  based  on  the  supposition  of 
hc  -v  H,  rather  plausible  core  configurations  seem  to  have  been 
drawn,  and  more  experiments  are  needed  to  verify  the  possible 
appearance  of  flow  subregimes  in  the  intermediate  flow  regime  of 
concern  here. 

(ii )  Pr  >>  1 

For  Pr  >>  1,  globally  the  core  flow  structure  was  shown  to  be 
parallel  in  both  Case  I  and  Case  III.  In  this  case,  due  to  the  high 
viscous  shear,  it  is  not  likely  that  the  development  occurs  in  such 
a  way  that  first  the  parallel  core  flow  structure  in  Case  I  becomes 
the  non-parallel  core  flow  structure  in  the  intermediate  flow  regime 
and  then  the  non-parallel  core  flow  structure  reverts  to  the  parallel 
core  flow  structure  in  Case  III.  Thus  the  parallel  core  flow  struc¬ 
ture  is  still  expected  in  the  intermediate  flow  regime  for  Pr  »  1. 
Actually  in  the  experiment  done  so  far  for  Pr  »  1  for  the  present 
configuration,  e.g.,  Kamotanl  et.  al.  [33]  and  Simpkins  and  Dudderar 
[38],  the  parallel  core  flow  pattern  was  observed  for  the  most 
ranges  of  the  parameter  RaH  unless,  as  mentioned  earlier,  secondary 
cells  appear  near  the  end  in  the  boundary-layer  flow  regime. 
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( 1  i  1 )  Pr  <  1 


In  this  case,  the  appropriate  force  balance  is  that  between 
the  buoyancy  and  inertia  forces.  We  first  consider  the  force 
balance  in  the  end  region.  From  (4.30),  this  balance  can  be 
represented  as 


_L 

4 


From  (4.35),  (4.36)  and  (4.59),  we  then  obtain 


(4.59) 


(_H  ~  ~Tfznr  (4-6°) 

In  order  to  determine  the  length  scales  hc  and  £  explicitly, 
one  more  force  balance  is  required.  Since  the  balance  between  the 
buoyancy  and  inertia  forces  in  the  core  yields  the  same  relation  as 
in  (4.59),  it  does  not  give  any  more  information.  Assuming  the 
possible  interaction  of  the  inertia  forces  between  the  end  and  core 
region,  one  additional  force  balance,  that  is,  the  balance  between 
the  inertia  in  the  end  and  the  inertia  in  the  core,  could  be  in¬ 
ferred.  From  Eq.  (3.40),  this  balance  can  be  represented  as 


and  this  gives 


or  from  (4.34) 


(4.61) 


(4.62) 
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1 

(PrRaH)1/3 


1 


(4.69) 


In  this  case,  patterns  of  multi -cel Is  may  be  possible  in  the  core. 
However,  in  order  that  PrRa^  ^  1,  Pr  should  be  of  order  0(10"4)  or 
less  than  that  for  the  common  ranges  of  the  parameters,  A  £  0.1  and 
RaH  £  10^  such  that  Ra^A2  >>  1.  Thus,  in  general,  that  does  not 
seem  to  be  a  realistic  situation. 

If  we  suppose  h£  ^  H,  from  (4.60)  we  obtain 

ic  *  (PrRaH)*  H  (4.70) 


The  characteristic  quantities  VR  and  in  (4.35)  and  (4.36)  can 
then  be  represented  as 


vR  %  a(PrRaH)* 


(4.71) 


and 


'V, 


H 

(PrRaH)* 


(4.72) 


These  characteristic  quantities  are  similar  to  those  of  the  boundary- 
layer  flow  regime  for  Pr  <  1.  This  also  implies  that  the  flow 
characteristics  in  this  intermediate  flow  regime  are  closer  to  those 
in  the  boundary  layer  flow  regime. 

Since  A(PrRaR)*  £  1  in  the  intermediate  flow  regmie  for  Pr  <  1, 
from  (4.70)  we  find 

*c  i 

-v  A(PrRaHr  £  1 


(4.73) 
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and  for  A(PrRaH)*  <  1 


(4.74) 


This  also  implies  the  existence  of  flow  subregimes  in  the  core  and 
the  length  scale  of  i  will  indicate  the  types  of  such  flow  sub- 

w 

regimes. 

When  A(PrRaH)*  <  1,  from  (4.70)  we  may  write 


_1 _ 

A(PrRaH)J 


(4.75) 


where  n  denotes  an  arbitrary  integer  approximating  the  ratio  of  y-. 

c 

Except  for  extremely  small  values  of  Pr,  from  (4.70)  we  also  find 


h  i 

f  -  - — r  <  1 

c  (PrRaH)i 

From  (4.75)  and  (4.76),  since 


h 

I 


c 

c 


% 


An 


(4.76) 


(4.77) 


we  thus  find 

1  <  n  <  i  (4.78) 

which  indicates  that  patterns  of  multi-cells  are  possible  in  the 

2 

core  for  A  «  1. 
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4  tr 

When  A(PrRaH)  ^  1,  from  (4.70),  ^  *  1.  In  this  situation, 
no  apparent  flow  subregimes  are  expected  and  the  core  flow  pattern 

is  not  readily  clear. 

4  lc 

When  A(PrRaH)  >  1,  —  >  1  and  thus  no  flow  subregimes  are 
expected.  This  is  the  case  of  boundary-layer  flow  regime.  It  is 
evident  from  the  global  analysis  that  no  such  flow  subregimes  are 
expected  in  the  boundary-layer  flow  regime  for  Pr  <  1. 

In  the  above,  although  there  do  not  exist  any  experimental  data 
available  for  direct  comparison,  the  core  configuration  drawn  from 
the  supposition  of  hc  ^  H  seems  to  be  a  rather  more  plausible  situ¬ 
ation  than  that  obtained  from  the  additional  force  balance  between 
the  inertia  forces  in  the  end  and  core  region. 

Utech  et.al.  [45]  conducted  an  experiment  with  molten  tin  of 
Pr  -v.  0(10  ).  Their  configuration  is,  although,  similar  to  that 
considered  here,  it  differs  significantly  in  that  the  upper  surface 
is  free  surface.  Experiment  shows  a  pattern  of  multi -cel Is  in  the 
core  superimposing  on  the  generally  steady  flow  along  the  bottom 
and  top  of  the  boat.  However,  this  is  supposedly  due  to  the  verti¬ 
cal  heat  transfer  through  the  open  boat  top. 

In  summary,  in  the  intermediate  flow  regime,  the  core  flow  pat¬ 
tern  is  not  parallel  and  flow  subregimes  such  as  secondary  cells  or 
patterns  of  multi -cel Is  seem  to  be  possible  both  for  Pr  ^  1  and 
Pr  <  1,  whereas  the  parallel  core  flow  pattern  Is  still  expected 
for  Pr  »  1. 
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CHAPTER  V 


SUMMARY  AND  CONCLUDING  REMARKS 

Consideration  has  been  given  to  the  prediction  of  core  flow 
pattern  in  a  low  aspect  ratio  rectangular  enclosures.  Globally 
core  flow  characteristics  were  studied  according  to  the  three  dif¬ 
ferent  cases.  With  the  information  obtained  from  the  global  anal¬ 
ysis,  the  core  flow  pattern  was  also  examined  in  detail  to  predict 
the  possible  flow  subregimes  such  as  secondary  cells. 

In  the  analysis  of  Case  I,  the  balance  was  made  between  buoy¬ 
ancy  and  viscous  forces  in  the  core.  The  analysis  based  on  this 

2 

force  balance  is  valid  within  the  parametric  ranges  of  Gr„A  £  1 
2 

and  RaHA  £  1.  In  this  flow  regime,  the  c jre  flow  is  parallel  and 
the  core  temperature  distribution  is  either  purely  linear  or  linear 
and  stratified. 

In  Case  II,  the  analysis  is  based  on  the  balance  between  buoy¬ 
ancy  and  inertia  forces  in  the  core.  The  analysis  in  this  case  is 

2  2  1 

limited  to  the  parametric  ranges  of  GruA  >>  1,  Pr  <<  — — 5  such 
that  PrRa^A4  «  l.  The  core  flow  is  not  parallel  in  this  flow 
regime  while  the  core  temperature  distribution  is  purely  linear. 

In  both  Cases  I  and  II,  the  flow  is  driven  in  the  core  and  no 
flow  subregimes  are  expected.  For  Ra^A  >>  1,  however,  the  driving 
force  was  shown  to  exist  in  the  end  region  and  this  situation  was 
considered  as  Case  III. 
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In  the  analysis  of  Case  III,  the  balance  was  made  between  con¬ 
vection  and  conduction  in  the  end  region.  In  addition,  the  balance 

was  made  in  the  end  region  between  buoyancy  and  viscous  forces  for 
Pr  *  1  (including  Pr  >>  1)  and  between  buoyancy  and  inertia  forces 
for  Pr  <  1  (this  always  includes  the  case  of  Pr  <<  1). 

In  Case  III,  under  the  conditions,  ARaH*  >  1  for  Pr  %  1 

(including  Pr  >>  1)  and  A(PrRaH)*  >  1  for  Pr  <  1,  there  exist  dis¬ 

tinct  horizontal  thermal  layers  adjacent  to  the  horizontal  boundar¬ 
ies  in  the  core  and  the  temperature  distribution  outside  the 
horizontal  thermal  layers,  i.e.,  in  the  mid-core,  is  stratified. 

The  core  flow  is  driven  by  the  thermal  boundary  layers  in  the  end 
region,  and  for  Pr  >>  1  the  core  flow  pattern  is  parallel  while 
the  fluid  motion  remains  nearly  stagnant  in  the  mid-core  for  Pr  ^  1 
(i.e.,  Pr  is  about  1  or  slightly  greater  than  1)  and  Pr  <  1.  In 
this  flow  regime  (called  the  boundary- layer  flow  regime),  flow 
subregimes  such  as  secondary  cells  are  expected  to  appear  near  the 
end  region  for  Pr  >>  1  under  the  condition  APr  %  RaH*,  while  in  the 

core  flow  subregimes  are  expected  for  Pr  ^  1  (except  Pr  £  1)  under 
1  Pr^ 

the  condition  - — ,  <  A  <  . 

Ran*  RaH 

When  ARa^*  ^  1  for  Pr  £  1  (including  Pr  »  1)  and  when 
A(PrRaH)*  *  1  for  Pr  <  1  ,  the  flow  regime  corresponds  to  the  so 
called  intermediate  flow  regime.  In  this  flow  regime,  the  core 
flow  pattern  is  not  parallel  and  flow  subregimes  such  as  secondary 
cells  or  patterns  of  multi -cel Is  seem  to  be  possible  in  the  core  for 
Pr  -v  1  and  Pr  <  1  whereas  a  parallel  flow  core  pattern  is  still 
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expected  for  Pr  »  1.  Summary  of  the  results  of  the  present 
analysis  is  given  in  Table  1. 

By  comparisons,  predictions  of  core  flow  pattern  made  herein 
show  good  agreement  with  the  existing  experimental  results  for 
Pr  -v  1  and  Pr  >>  1,  although  the  data  are  not  extensive.  For  Pr  <  1, 
no  experimental  data  for  the  present  configuration  are  available, 
as  yet,  for  direct  comparison.  More  experiments  are  needed  to 
verify  the  present  predictions  in  order  to  understand  the  physics  of 
the  core  flow  pattern  more  clearly  and  to  indicate  the  validity  of 
the  type  of  analysis  presented  herein.  Clearly  the  analysis  can  be 
generalised  for  other  configurations. 

Based  on  the  multiple  scales  technique,  global  prediction  of 
the  core  flow  pattern  is  generally  satisfactory.  Further  work  will 
be  required  to  define  the  flow  subregimes  better.  Nevertheless,  it 
is  thought  that  employing  the  ideas  of  multiple  scales  technique 
attempted  herein  can  surely  be  applied  to  resolve  many  other  complex 
probl ems . 
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APPENDIX  A 


Examination  of  the  situations  for  which  either  Eg.  (3.22)  or  Eg. 
(3.23)  could  be  satisfied. 


°  l3-22> 

O 

(3-23) 


(a)  »  «  0 

This  implies  no  flow  in  the  core  and  inappropriate  to  the 


current  situation. 

(b)  e  =  eU) 

In  this  case,  from  Eq.  (3.15) 

aLaj.  =  o 

->5 

(A.l) 

and  thus,  ^  =  ) 

(A. 2) 

Then,  from  Eqs.  (3.7)  and  (3.8) 

(A. 3) 

From  (A. 3),  we  find 

0  -  Const. 

(A.4) 

This  Implies  no  driving  force  in  the  core,  thus  it  is  inappropriate 
to  the  current  situation. 
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(c)  H  "  °*  iuy  *  0 


In  this  case,  ^  and  e  can  be  represented  as 


f  -  ft}) 

9  =  fi%)  t 


(A. 5) 


Substitution  of  (A. 5)  and  (A. 6)  into  Eqs.  (3.7)  and  (3.8)  yields 


f'<s>  *  f'a> 


(A. 7) 


Comparing  both  sides,  the  equality  can  be  made  when  each  side 
has  constant  value.  Denoting  the  constant  by  A^,  from  (A. 5),  (A. 6) 
and  (A.7),  we  find 


0  s  M  +  }<}> 

f  *  a, 


(A. 9) 


Eqs.  (A. 8)  and  (A. 9)  show  that  the  temperature  distribution  in 
the  core  is  linear  and  stratified,  while  the  core  flow  pattern  is 
parallel . 


(d)  i  -  6  (y) 


Considering  the  boundary  conditions  (2.9b),  we  expect  the 
functional  form  of 


f  =■ 


(A. 10) 
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and 

IjJ-  «  f'mktp 

(A. 11) 

or 

&  *  f  t  fr*>,  -  j 

(A. 12) 

such  that 

(A. 13) 

Substituting 

(A. 10)  and  (A. 12)  into  Eq.  (3.7),  and  from  Eq. 

(3.8) 

we  have 

(A. 14) 

we  separate  the  variables  in  (A. 14)  as 

-  JJ1L.  =  J —  f  -  1&-  ] 

PW  »■  ®  J 

This  equation  is  separable  if 

$jf  -  «*«•  -  ** 


or 


Pit)  -  const. 


(A. 15) 


(A. 16) 
(A. 17) 


When  p(y)  =  const.,  from  (A. 12) 
0  *  £<■$> 


(A. 18) 


This  situation  was  previously  treated  in  case  (b)  and  it  was  shown 
to  be  inappropriate  to  the  current  situation. 

For  the  case  (A. 16),  denoting  the  suitable  separation  constant 
by  k,  from  (A. 15)  and  (A. 16) 


-  kptl)  + 


(A. 19) 


105 


Let  us  look  over  (A. 19)  according  to  the  values  of  k  and 
(i )  k.  A?  arbitrary  and  non-zero 

Considering  the  symmetry  property  of  in  (2.10a), 


fit.*)  *  /-*> 

(A. 20) 

From  (A. 10)  and  (A. 20), 

(A. 21) 

We  thus  find 

(A. 22) 

and 

rs 

ar\* 

\ 

1) 

o\> 

(A. 23) 

Also  considering  the  symmetry  property  of  6  in  (2.10b), 

*)  -  / -  »-*> 

(A. 24) 

From  (A. 12)  and  (A. 24),  we  then  have 

t  flS)  ~  /-  f'c-UPd'*)- t(i-l) 

(A. 25) 

Since  f-1(c)  *  f1  (1-c)  from  (A. 22),  we  find 

p<*>  *  - 

(A. 26) 

and 

*  '  '  ft/- J> 

(A. 27) 

From  (A. 23) 

(A. 28) 

From  (A. 19)  and  (A. 28),  we  thus  have 


(A. 29) 
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Since  from  (A. 26)  and  (A. 29) 

PP<t)  -  -  K  pit) 

we  find 

PCD  =  0 

Then,  from  (A. 12) 

0  -  fill) 


(A. 30) 


(A. 31) 


(A. 32) 


This  is  identical  to  the  case  (b)  which  was  shown  to  be  inappropri¬ 
ate  to  the  current  situation 


(ii)  A2  =  0 


In  this  case,  from  (A. 16) 

12a  =  o 

tt%) 

Since  e  is  finite,  from  (A. 12),  f(c)  t  0.  Thus 

r«>  -  o 


and 


fit)  =  const. 


From  (A. 19) 


fci)  =  UPC}) 


(A. 33) 

(A. 34) 
(A. 35) 

(A. 36) 


Considering  the  symmetry  properties  of  e  and  t  from  (A. 26)  and  (A. 28) 
we  find 


PH)  -  0 


(A. 37) 


r 

i 

i 
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Then,  from  (A. 12)  and  (A. 33)  e  =  const.,  and  this  is  inappropriate 
herein. 


(iii)  k  °  0 


In  this  case,  from  (A. 19),  g1v(y)  =  A^,  and  from  (A. 15)  and 


(A. 16) 


JjlL  =  0 


(A. 38) 


Since  f(c)  i  0  for  finite  e  in  (A. 12), 

'  0  (A. 39) 

and  f(f)  -  const.  =  A3  (A. 40) 

Then,  from  (A. 16)  and  (A. 40) 

*  Ax  fd)  -  Ax  A,  =  (A.  41) 

and  £(*1  *  A.%  +  A,  (A.42) 

where  A&  is  an  arbitrary  constant. 

Substituting  (A. 40)  and  (A.42)  into  (A. 10)  and  (A. 12),  we  finally 
have 

t  s  (A. 43) 

and  6  =  p( V  +  +  A,  (A. 44) 

These  are  the  same  results  obtained  in  the  analysis  of  case  (c). 

(e)  ||-  dy  ■  0 


r 

l 


From  the  symmetry  properties  of  y  and  e. 
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f(f)  =  fu-v 


m  -  il\ 


(A. 45) 


(A. 46) 


Since  both  ij;(y)  and  —  have  the  same  sign  between  y  =  0  and  1, 

ay 

there  exists  no  possibility  of  the  value  of  the  integral  dy 

3  A 

to  be  zero,  unless  t  =  0  or  =  0  (i.e.,  e  =  5(c))  in  which  case 

ay 

it  was  already  shown  to  be  inappropriate  in  the  case  (a)  and 


case  (b),  respectively. 


/f\  it  ae  _  ,  a  6 

1 f 1  3;  ay  ~  3s  3y 


(A. 47) 


Considering  the  boundary  conditions  (2.9b),  \ and  e  can  be 


expressed  as 


f  -  fd) 

9  = 


(A. 48) 
(A. 49) 


Substituting  (A. 48)  and  (A. 49)  into  (A. 47)  and  dividing  each  side 
by  fU)g(y)F(t)G(y),  we  then  have 


fv i>  _  _  p*») 

f  ~  F<S> 
Integration  of  (A. 50)  yields 

'inf  *  ~  in F  +  ^n.A4 


(A. 50) 


(A. 51) 
(A. 52) 


where  Ag  is  an  Integration  constant. 
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From  (A. 49)  and  (A. 52),  we  find 

0  *  A‘  ^7  (A. 53) 

From  (A. 24)  and  (A. 53) 

(A-54> 

Since  from  (A. 20)  and  (A. 48) 

f(5>  -  (A. 55) 

from  (A. 54)  and  (A. 55),  we  then  have 

<*■«> 

or  +  4(/-fr>)  =  ff5)  (A. 57) 


Comparing  both  sides,  the  equality  can  be  made  only  when  each  side 
has  a  constant  value.  Denoting  the  constant  by  A7,  from  (A. 48)  and 


(A. 57)  we  find 

f  *  Aj  fr(jf)  (A. 58) 

From  (A. 47)  and  (A. 49),  this  gives 

F'(1)  *  o  (A. 59) 

and  thus  Fd)  -  const.  (A. 60) 

or  Cnf)  =  const.  (A. 61) 


The  situation  (A. 60)  represents  no  driving  force  in  the  core  and  the 
situation  (A. 61)  reduces  to  Case  (b).  This  case,  therefore,  does  not 
show  any  possible  core  configurations. 


APPENDIX  B 


Examination  of  the  situations  for  which  either  Eg.  (3.78)  or  Eg. 
(3.79)  in  Case  II  could  be  satisfied. 


or 


(a)  &  =  0 


(3.75) 

(3.76) 


This  means  no  flow  in  the  core  and  it  is  inappropriate  to  the 
current  situation 


(b)  6  =  eU) 


In  this  situation,  from  Eq.  (3.74) 


and  thus. 


22-ZL  -  0 
f  * 


Then  from  Eqs.  (3.67)  and  (3.69),  for  tc  <<  1 

M.  =  * 
dS 


and 


0  =  const. 


(B.  1 ) 
(B.2) 


(B.3) 

(B.4) 


This  implies  no  driving  force  in  the  core,  thus  it  is  inappropriate 
to  the  current  situation. 


Ill 


In  this  situation,  <i>  and  e  can  be  represented  as 

f  -  fd)  (B. 5) 

and  0  =  fw)  -t  $<»)  (B.6) 

Substitution  of  (B.5)  and  (B.6)  into  Eq.  (3.69)  yields 

0  -  (B.7) 

Thus,  const.  =•  f(%)  (B.8) 

Then,  from  (B.6)  and  (B.B) 

9  *  const,  t  (B.9) 


This  implies  the  stratification  cf  core  temperature  distribution 
and  no  driving  force  in  the  core,  thus  it  is  inappropriate  to  the 
current  situation. 

(d)  !  ‘  a7  =  F(y) 


In  this  situation,  as  in  the  analysis  of  case  (d)  in  Appendix 
A,  if  we  consider  the  boundary  conditions  (2.9b),  e  and  c  can  be 
expressed  as 


t  - 

(B. 10) 

and 

ff-  *  fU)h(t) 

(B.ll) 

or 

p  *  +  ]*<*>*} 

0 

(B. 12) 

so  that 

f| $■  =  ftfihCf)  *  /?(*) 

(B. 13) 
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Substituting  (8.10)  and  (B.12)  into  Eq.  (3.69)  and  from  Eq.  (3.67), 

2 

for  A  «  1,  we  have 

ftlitfDlt'tvfw-itvFii)}  *  +  (B.  14) 


After  arranging. 

jlf; )  - 

f ''12-  1  -  i 

(B.15) 

This  equation  is  separable  if 

pit)  =  const. 

(B. 16) 

°r  iMtk  = 

(B.17) 

When  p (y )  =  const.,  from  (B.12)  e  =  e(;).  This  situation  is 
shown  to  be  inappropriate  to  the  current  situation  in  case  (b). 

For  the  case  (B.17),  denoting  the  separation  constant  by  k  in  (B.15), 


from  the  second  expression 

f(  57  -  K  =  const.  (B.  18) 

Then  from  (B.10)  and  (B.12), 

f  -  (B.  19) 

and  e  =  K*  P<f)  t  f(S>  (B. 20) 

Substitution  of  (B.18)  into  (B.14)  yields 

j-U)  *  const.  (B.21) 


thus. 


f(1)  =•  const. 


(B.22) 


This  is  also  consistant  with  the  result  which  can  be  obtained  from 
(B.17)  and  (8.18).  From  (B.20)  and  (B.21),  e  then  becomes 

9  *  K  ?<*>  +  const.  (B- 23) 

This  shows  that  the  core  temperature  distribution  is  stratified  and 
it  does  not  provide  any  driving  force  for  the  core  flow.  It  is  thus 
inappropriate  to  the  current  situation. 


By  the  same  reasoning  of  case  (e)  in  Appendix  A,  this  does  not 
show  any  possible  core  configurations. 


By  the  same  procedure  of  case  (f)  in  Appendix  A,  this  case 
reduces  to  the  case  (c)  which  does  not  show  any  possible  core 
configurations  appropriate  to  the  current  situation. 


APPENDIX  C 


Consideration  on  the  Possible  Core  Configurations  from  Egs.  (3.110) 


From  Eq.  (3.110),  outside  the  horizontal  thermal  layer  two 
possible  core  temperature  profiles  can  be  obtained,  i.e.,  ec  = 
const,  and  e  *  e  (y).  Amont  them  e  *  e  (y)  is  a  possible  one, 

v  L  v  V, 

thus  we  are  going  to  examine  whether  ec  *  const,  is  also  possible 
core  temperature  profile.  For  the  analysis,  we  develop  simple 
integral  relation  which  can  describe  the  core  configurations. 

By  introducing  the  transformation  (C.l) 


/ 


(C.l) 


into  Eq.  (3.100),  the  core  energy  equation  within  the  horizontal 
thp'-mal  layer  becomes 


<‘-2> 

Rewriting  Eq.  (C.2),  we  have 

'  A*H~)  ~  (C-3) 

Integration  of  Eq.  (C.3)  with  respect  to  y  between  0  and  i.e., 
between  the  horizontal  wall  and  mid-core  region  (i.e.,  outside  the 
horizontal  thermal  layer),  yields 

‘  JfM f (c-4) 
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The  left  side  is  Integrated  by  parts  and  we  then  have 


and  after  arranging 


where 


(C.5) 


(C.6) 


«  r\„m  ,  t  =  1 19.. 

and  subscripts  c  and  w  represent  the  horizontal  wall  and  mid-core, 
respectively. 

By  applying  the  boundary  conditions  (2.9b)  at  horizontal 
boundaries,  Eq.  (C.6)  reduces  to 


-IjM" (c-e) 

Carrying  out  the  derivatives  on  left  side  and  after  rearranging, 
we  finally  have 


tff  -  &(YIH  -  3f-|?„+  (c.») 

The  integral  equation  (C.9)  describe  the  core  configuration. 
Substitution  of  e  *  const,  into  Eq.  (C.9)  then  yields 

t 
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For 


<<  1,  Eq. 


(C.IO)  reduces  to 


(C.10) 


Irl'H m  * 0 


(C. 11) 


Let  us  examine  what  the  situations  for  which  either  Eq.  (C.ll) 
or  Eq.  (C.12)  could  be  satisfied  would  imply. 

(a)  0 

This  implies  no  flow  within  the  horizontal  thermal  layer  and 
thus,  no  flow  in  the  core.  It  is  thus  inappropriate  to  the  current 
situation  because  there  is  no  heat  transfer  between  the  two  ends. 


(b)  e  =  const. 


Since  0„  =  const,  too,  the  entire  core  remains  at  constant 
c 

temperature.  This  means  that  no  heat  is  transferred  across  the  core 
by  convection,  thus  it  is  inappropriate  to  the  current  situation. 


(C)  .  0.  Me  ■  o 

3;  3c3y 


In  this  case  ^  and  e  can  be  represented  as 
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t  *  tC*>  (C.  13) 

and  0  *  ♦  }(})  (C.  14) 

when  y  -  «  ,  from  (C.14) 

el  s  6C  =  «■  *<«>)  *  . onst .  (C.  15) 

Since  &c  is  considered  to  be  const,  herein,  this  is  not  appropriate 
to  the  current  situation. 


Considering  the  boundary  conditions  (2.9b),  we  expect  the 
functional  form  of 


(C. 16) 

and 

If"  =  f"(  5)K(>) 

(C. 17 ) 

so  that 

s  =  F(*> 

(C. 18) 

In  this 

situation,  from  (C. 17) 

0  =  t  frC*) 

(C. 19) 

with 

(C.20) 

when  y  ■*  »  ,  from  (C .  19) 
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0|?„„  *  6C  »  t  *  const.  (C.21 ) 

Thus,  this  is  also  inappropriate  to  the  current  situation. 

(e)  fo  *  4  dy  =  0  (C.22) 

_ 33g _ 

Considering  the  symmetry  properties  of  *  and  o,  as  can  be 

inferred  from  the  case  (e)  of  Appendix  A,  c  and  —  are  either 

3y 

increasing  or  decreasing  between  the  horizontal  wall  and  the  mid¬ 
core  region.  Thus,  there  is  no  reason  that  the  integral  in  (C.22) 
becomes  zero,  unless  ■„  =  0  or  ~  =  0,  i.e.,  -  =  ?(;)  or  const., 
in  which  case  it  was  shown  to  be  inappropriate  to  the  current  situa¬ 
tion  in  case  (a)  and  case  (b),  respectively. 

From  the  considerations  made  above,  it  is  shown  that  e  = 

const,  is  not  a  possible  core  temperature  profile  herein.  However, 

2  2 
from  Eq.  (C.10),  when  A  is  not  negligibly  small,  i.e.,  for  A 

0(1),  the  right  expression  in  Eq.  (C.10)  is  not  negligible  and  Eq. 

(C.10)  itself  can  h*  satisfied  for  a  suitable  distribution  of  f 

within  the  horizontal  thermal  layer.  Thus,  ec  =  const,  in  this 

situation  can  also  be  a  possible  core  temperature  profile.  This 

means  that  heat  is  transferred  by  convection  across  the  cavity 

through  the  horizontal  thermal  layer  while  the  mid-core  region 

remains  at  constant  temperature. 


APPENDIX  D 


Modified  Characteristic  Stream  Function 

In  the  boundary  layer  equations,  in  normalizing  the  velocity 
components  within  the  thermal  or  flow  boundary  layers  we  use  the 
same  characteristic  velocity,  UR,  because  the  characteristic 
velocity  is  of  the  same  order  of  magnitude  within  both  boundary 
layers.  On  the  contrary,  in  vorticity  transport  equations  of 
the  boundary  layer,  since  the  characteristic  stream  function  is 
represented  as 


where  t  is  the  characteristic  length  ,  even  for  the  same  value 
of  UR,  the  characteristic  stream  function  changes  according  to  the 
characteristic  length  scale,  1,  and  thus  one  should  be  careful 
in  normalizing  the  vorticity  equations.  Actually,  in  the  case  of 
large  Prandtl  numbers,  the  flow  boundary  layer  extends  beyond  the 
thermal  boundary  layer,  thus  the  characteristic  stream  function 
differs  within  each  layer  and  it  needs  to  be  determined  according 
to  the  corresponding  boundary  layers.  For  low  Prandtl  numbers, 
however,  although  the  thermal  boundary  layer  extends  beyond  the  flow 
boundary  layer,  since  the  main  body  of  fluid  flows  within  the 
flow  boundary  layer,  the  characteristic  stream  function  will  be  of 
the  same  order  within  both  boundary  layers.  In  Eq.  (3.89  ),  the 


119 


120 


estimate  of  the  characteristic  stream  function,  yr,  was  made  for 
the  thermal  boundary  layer,  say  6^,  for  Pr  K  1  (including  Pr  >>  1). 
Based  on  the  above  arguement,  for  the  same  Prandtl  numbers,  we  now 
determine  the  characteristic  stream  function,  ?R,  for  the  flow 
boundary  layer. 

From  (D.l),  UR  can  be  written  as 


(D.2) 


where  ?R  =  —  >fR,  is  the  characteristic  stream  function  modified 
for  the  flow  boundary  layer,  6^.  To  estimate  VR,  we  replace  vR 
by  Tr  in  Eq.  (2.13)  and  introduce  the  multiple  scales  for  the 
flow  boundary  layer,  6^,  as 

C  =  X,  n  -  (D.3) 

x 

6f 

where  j—  (D.4) 


Then  substituting  the  derivatives  of  (0.3)  into  Eq.  (2.13)  and 
balancing  the  inertia  and  viscous  forces  within  6^  in  the  end,  we 
have 


¥i7 


(D.5) 


and  from  (0.2)  and  (0.5),  we  find 


(D.6) 


For  since  from  (3.85) 


j 
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. - . H 

* 

e  ^  SA 

X  vR 

(D.7) 

t 

i 

1  • 

i  j 

from  (D.6)  and  (D.7),  we  have 

1: 

6t 

rx  '  Pr  67  cx 

(D.8) 

H 

i 

or  —  -v.  Pr 

ex  6f 

(D.9) 

\ 

i 

i 

! 

Then,  from  (3.52),  (D.4)  and  (D.9) 

J 

rv  6.  6. 

1 

T“  ^  ^  Pr  t — 

£x  6f  6f 

(D. 10) 

I 

i 

®f  i 

and  thus,  — -  ^  Pr1 

(D.ll) 

I 

i 

6t 

l 

f 

1. 

Since  from  (3.91) 

1 

i: 

6  ^  H 

1  Ra  1 

KaH 

(D.12) 

1  ’ 

I  ' 

i 

from  (D.ll)  and  (0.12),  we  find 

) 

| 

i. 

r 

6f  -v  H 

(D.13) 

i 

j, 

i. 

_  i  . 

and  thus,  vR  =  vR  <v  a  Pr*  RaH* 

(D. 14) 

1  J 

k  1 

l 

With  this  modified  characteristic  stream  function,  ?R 

,  in  (D. 14) , 

!  : 
1  ■ 

j 

[ 

C 

i 

] 

the  estimate  of  the  horizontal  viscous  layer  was  made 

in  (3.119) t 

1 

j 

j 

j 

1 

t 

i 

] 

• 

i 

a 

i 

i 

’i 
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